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ABSTRACT 


This report details the development and provides the 
documentation for custom computer software applications that 
evaluate antenna parameters. The applications are written in 
Mathcad 3.1 for the following antenna types: linear, planar, 
and circular arrays; folded dipole; caged dipole; parabolic 
reflectors with helical and spiral feeds. Inputs to the 
applications are limited to the antenna’s physical dimensions. 
In some cases, ground parameters are required. 

The chapters are structured to provide the user with the 
necessary information needed to use and interpret the software 
for each antenna type. The software applications are provided 
aS appendices and give examples of each antenna type. 

Outputs of the applications provide various numerical and 
performance predictions, as well as far-field radiation 
patterns. The results computed are consistent with 
predictions provided in applicable publications, as well as 
those calculated by numerical antenna analysis programs such 


as ELNEC. 
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I. BACKGROUND AND PURPOSE 


This thesis details the development of user friendly 
mathematical applications capable of computing the radiation 
pattern and other pertinent antenna parameters of an antenna 
or antenna system based on available information, primarily 
dimensional information obtained from photographs. The 
applications developed are compatible with any IBM personnel 
computer using MS-DOS version 3.2 or higher and a math 
coprocessor. The applications are written in Mathcad version 
3.1 engineering software and can be run using either Mathcad 
3.1 or Mathcad version 4.0. These applications are intended 
for use by the Naval Maritime Intelligence Center 
(NAVMARINTCEN) . 

As already mentioned, user inputs are limited to antenna 
dimensions based on physical measurements primarily obtained 
from photographs, although other source data can be used. The 
applications are flexible to the extent that other parameters 
may be estimated and used as inputs to increase the accuracy 
of the computations to gain better insight into the antenna’s 
performance. The Mathcad applications provide various 
performance predictions as well as a graphical representation 
of the antenna’s far-field radiation pattern. The necessary 


background information needed to interpret the application’s 


formulas and displays are provided in the corresponding thesis 
chapters. 

Dietrich [Ref. 1] and Gerry [Ref. 2] have completed the 
first and second reports of this project, respectively. This 
thesis is the last in a series of three reports intended to 


fulfill the NAVMARINTCEN statement of work. 


II. INTRODUCTION 


The increasingly sophisticated design by foreign 
countries of antennas warrants the development of improved 
analytical methods to provide timely and accurate technical 
assessment of these antennas. Through the process of reverse 
engineering, foreign communication, navigation, 
Identification, Friend or Foe (IFF), and radar antennas are 
evaluated to determine their capabilities, limitations, and 
vulnerabilities. Having limited information available 
concerning the parameters of an antenna, usually restricted to 
photographs that provide only the type and physical dimensions 
of the antenna, current methods of antenna evaluation are slow 
and tedious. 

With the availability of high-speed personal computers 
matched with sophisticated off-the-shelf engineering software, 
the assessment of antennas using only photographic information 
may now be obtained rapidly. The objective of this report is 
to provide NAVMARINTCEN with user friendly Mathcad software to 
achieve this type of performance in their evaluation of 
antenna systems. 

As with any computer analysis of antennas, the accuracy 
of the results depend upon the sophistication of software 
used, the complexity of calculations, the reliability of input 


data obtained from photographs, and the accuracy of any 


required input data that must be estimated. This report and 
the Mathcad applications are written with existing engineering 
equations for the type of antenna analyzed and inputs are 
limited to physical dimensions. Additional information such 
as ground effects and characteristic impedance may be included 
to increase the accuracy of the analysis. To this end, this 
report should provide an initial understanding as to the 
capabilities of the antenna system in question. 

The chapters of this report document the equations and 
assumptions used for each specific antenna type and is written 
as a comprehensive reference for the software. Copies of each 
application are included as appendices to provide the user 


with a printed illustration of the software. 


III. THE UNIFORM LINEAR ARRAY ANTENNA 


Arrays are highly directional antenna system with narrow, 
steerable beams and low side lobes used for long distance 
communication and radar systems. The linear array is one 
geometrical configuration that consists of radiating elements 
lying along a straight line. As a result of the dependance of 
inter-element spacing (d) to the wavelength (A), linear arrays 
are inherently narrowband antenna systems. 

In this Mathcad application uniformly spaced, equally 
excited, isotropic point sources aligned on the z-axis are 
used as the radiating elements. In using isotropic point 
sources, polarization of the array is not calculated and an 
assumption is made that the antenna input resistance (R,,) is 
equal to the radiation resistance (R,). This assumption 


results in calculations for gain (G) that are idealized 


values. The following inputs are needed to implement the 
software: 

N = number of isotropic radiating elements 

d = inter-element spacing between adjacent elements 

6, = direction of main lobe : 

f = frequency of interest “ 

I. = antenna feed current ~ 


oO 


Z 


oO 


. . . * 
characteristic feed impedance 


The first two inputs are obtained from visual data such as 


photographs, and the ~ indicates input parameters that are 
either known or are estimated. The linear array geometry is 


illustrated in Figure 3.1. 


FIGURE 3.1 Geometry of N-element linear array of 
isotropic point sources 


The total electric field of a linear array in the far- 
field, neglecting mutual coupling between adjacent elements, 
is determined by the vector addition of the fields of 
individual radiating elements. It can be shown that the total 
field is equivalent to the product of the field of a single 
element, selected at a reference point, and a factor which is 


referred to as the array factor (AF): 


E(total) = (E(element)] x [AF] (3.1) 


In (3.1), E is the vector electric field intensity, while AF 
is a scalar quantity. This result is referred to as pattern 
multiplication for arrays of identical elements. The array 
factor is a function of the geometry of the array and the 
element excitation amplitude and _ phase. For uniform 
excitation, by varying the inter-element spacing (da) and/or 
the phase (8) between the elements, the characteristics of the 
array factor, and consequently the total electric field, of 
the array can be controlled. [Ref. 3: pp. 204-207] 

This Mathcad application is written for an N-element 
linear array with equally spaced, uniformly excited isotropic 
radiating elements. In addition, element phase is assumed to 
vary linearly along the length of the array. The radiation 
pattern of a specific element type is neglected in (3.1) since 
in normal usage it will have little effect when the array 


consists of a large number of elements. However, the user of 


this software should be able to use information regarding the 
type of element and the orientation of the elements to narrow 
the range of possible 6,’s. The normalized array factor (AF) 


is (Ref. 3: pp. 212-214]: 


sin( yw (8)| 
AF, (0) = + (dimensionless) (3-2) 
N sin(>¥ (6) 
2 
where 
w (6) = kdcos(68) + B (radians) (3.3) 
B = -kdcos(8,) (radians) (3.4) 
In (3.4), @ is the direction of the maximum array factor 


° 

value (i.e., for a broadside array 0, = 7/2, for an end-fire 
array with the main beam directed at 180°, 6, = 7, and for a 
phased array with the main beam directed to +60°, 8, = 7/3), 
B is the phase excitation difference between the elements, k 
is the wavenumber (27/4), and ad is the inter-element spacing. 

With the linear array aligned along the z-axis, the array 
factor is a function of 8 but not of ¢. The visible region of 
the array in terms of 6, (8), and the inter-element spacing 


(a/A) is: 


0<O<n (radians) (3.5) 


B-kd < w(8) < Btkd (radians) (3.6) 


d/A determines how much of the array factor appears in the 
visible region as defined by (3.5) and (3.6). Although d/A 
provides guidance in determining the frequency to select for 
this Mathcad application, it should be noted that most arrays 
are designed with the inter-element spacing less than one 
wavelength, usually close to a half-wavelength. Thus, linear 
arrays are characteristic narrowband antenna systems. ([Ref. 
4: p. 123] 

An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r21.6A (meters) (3.7) 
r 25 (Na) (meters) (3.8) 
rs 2a (meters) (3.9) 


The array length, Nd, in (3.8) and (3.9) is the maximum 
dimension of the antenna (the array length is assumed to 
include a distance of d/2 beyond each end element [Ref. 5: 
p.55]). The minimum distance to the far-field is found by 
selecting the largest value of (3.7), (3.8), and (3.9). 
Radiation intensity (U(8)) is related to the power 


radiated in a given direction and is independent of the 


distance to the observation point [Ref. 3: p. 27; Ref. 4: p. 
33]. For linear arrays U(8) is given by [Ref. 3: pp. 229- 
233]: 


U(08) = (AF,)? (w/ solid angle) (3.10) 


Radiated power (P..,) is the total power radiated by the 
antenna and is obtained by integrating U(8#) over a surface 
surrounding the antenna [Ref. 4: p. 33]. Since the geometry 
of the linear array is aligned on the z-axis and the array. 


factor is independent of ¢, P.., is thus [Ref. 3: pp. 229-233]: 


Prag = 27 [790 8) sin (8) a0 (W) (3.11) 


Directivity (D,) is the maximum value of directive gain 
where directive gain is the ratio of U(@) in a specific 
direction to the average radiation intensity [Ref. 4: pp. 34- 
36]. Hence, D. is given by [Ref. 3: pp. 229-233]: 


fe) 


4u U. ; . 
D, = ——™ (dimensionless) (3.12) 


In (3.12), U,,, equals unity and occurs at 8, for the uniform 
linear array. 

Effective isotropic radiated power (EIRP) is a commonly 
used communications term that is defined as the product of 


antenna gain and total power accepted by the antenna from the 


transmitter. EIRP is determined as (Ref. 4: p. 62]: 


10 


EIRP = PragD,  (W) (3.13) 


Assuming that the power delivered to the array consisting 
of N-isotropic point sources is equivalent to the radiated 
power, we can estimate the radiation resistance (R,.) of the 


antenna as [Ref. 3: p. 55]: 


2 (Prag) 
= (Q) 3.14 
iz, |? ae 


Ir 


Using the same logic, we obtain the input resistance (R,;,) of 
the array made up of isotropic point sources as equivalent to 
the radiation resistance: 


(Q) (3.15) 


The exact gain (G) of the array cannot be calculated with 
the data assumed, but an idealized value can be computed based 
on the assumption made to determine R,, and the assumptions 
made to calculate total antenna efficiency (€,). The total 


antenna efficiency is given by [Ref. 3: pp. 44-45]: 


€é, = €,€,,=1-|P|? (dimensionless) (3.16) 

where 
E.g = 1 (dimensionless) (3.17) 
e,=1-|[|? (dimensionless) (3.18) 
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[ = —in “eo (dimensionless) (3.19) 


In (3.17), € is the conduction and dielectric efficiency 
which is unity for an ideal lossless antenna. In (3.19), Z, 
is the characteristic impedance of the transmission line and 
I is voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(3.18) and the ideal gain for a linear array is [Ref. 3: pp. 
43-44): 


(dimensionless) (3.20) 


G(dB) =10log,,(e€,D,) (dB) (3.21) 


As a result of constructing the linear array with 
isotropic point sources, the polarization loss factor (PLF) 
and antenna polarization are not calculated. Therefore, the 
maximum effective aperture (A,,) of the array when used as a 
receiver is calculated assuming negligible polarization 


mismatches. The A, is then [Ref. 3: p. 63]: 


a, = S14) per (m3) (3.22) 
4m 


where 


PLF =1 (dimensionless) (3.23) 
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The assumptions used in finding R, and A, result in an 
idealized value in calculating the effective height (h,). The 


effective height is [Ref. 6: p. 42]}:] 


RA 
Pio 2 | (meters) (3.24) 


where n, is the characteristic impedance of free space (1207). 

Two displays are produced in this Mathcad application. 
The first display is the rectangular representation of the 
magnitude of the array factor from -7 to 7. The second 
display is the polar representation of the magnitude of the 
array factor. 

The Hansen-Woodyard end-fire array was not addressed in 
this application nor were arrays of uniform spacing but 
nonuniform amplitude (i.e., Dolph-Tschebyscheff array). In 
addition, nonuniformly spaced, uniformly excited arrays were 
not addressed by this application. 

The results obtained with this application are identical 
to those found in [Ref. 3: pp. 216-240} with respect to array 
factor patterns and directivity. Table 3.1 is a comparison of 
data from [Ref. 3} for a 10-element equally excited, uniformly 
spaced, end-fire linear array while Figure 3.2 is a comparison 
of the array factor pattern of the same array. Of note, the 
array factor patterns from this application and the patterns 


from [Ref. 3] are identical as expected. 


i3 


Table 3.1 End-fire Array Data Comparison 


ANTENNA REFERENCE CALCULATED 


DIMENSIONS DIRECTIVITY (D,) DIRECTIVITY (D,) 
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FIGURE 3.2 Array factor pattern of 10-element equally 
excited, uniformly spaced end-fire linear array with d = 
A/4 m, A. = 0°, and B = -kd 


(9) 
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Table 3.2 is a comparison of data from (Ref. 3] fora 10- 
element equally excited, uniformly spaced, broadside linear 
array, and Figure 3.3 is a comparison of the array factor 
pattern of the same array. The array factor patterns from 


this application and the patterns from (Ref. 3] are identical. 


Table 3.2 Broadside Array Data Comparison 


ANTENNA REFERENCE CALCULATED 


DIMENSIONS DIRECTIVITY (D,) DIRECTIVITY (D,) 
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FIGURE 3.3 Array factor pattern of 10-element equally 
excited, uniformly spaced broadside linear array with d@ = 
4/4 m, 8, = 7/2, and B = 0 
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IV. THE PLANAR ARRAY ANTENNA 


Planar arrays, frequently used long distance 
communications and radar systems, exhibit characteristics 
analogous to those of linear arrays but have additional 
capability to control and shape the radiation pattern. Planar 
arrays are more versatile than linear arrays, providing 
symmetrical patterns with lower side lobes, and can be used to 
scan the main beam of the antenna array toward any point in 
space [Ref. 3: p. 261]. As with linear arrays, planar arrays 
are narrowband antenna systems due to the dependance of the 
inter-element spacing (d, and d,) on the wavelength (A). The 
following inputs are required analyze a planar array: 

M = number of isotropic radiating elements in the 

x-direction 

N = number of isotropic radiating elements in the 

y-direction 

d, = inter-element spacing between adjacent elements in 

the x-direction 
dad, = inter-element spacing between adjacent elements in 

the y-direction 
= frequency of interest * 


* 


direction of main lobe at 8 = oe 


f 
6, 


ce) 


Po direction of main lobe at ¢ = @ 
I 


antenna feed current “* 
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Z, = characteristic feed impedance e 


The first four inputs are data obtained primarily from visual 
data such as photographs, and the * indicates input parameters 


that are either known or estimated. The planar array geometry 


is illustrated in Figure 4.1. 


FIGURE 4.1 Geometry of M x N element planar array 
of isotropic point sources 


19 


This Mathcad application is written for an M x N planar 
array with uniformly spaced, equally excited, identical 
radiating elements positioned on the x-y plane. This 
application may also be used to approximate the array factor 
and field patterns for arrays that have elliptical or circular 
geometries by assuming an ellipse is equivalent to a rectangle 
or a circle is equivalent to a square, respectively [Ref. 6: 
p.187]. Since isotropic point sources are assumed as array 
elements, the polarization of the array is not calculated, and 
an assumption is made that the antenna input resistance (R,,) 
is equivalent to the radiation resistance (R,). This 
assumption results in values for gain (G) that are idealized. 
As with the linear array, the radiation pattern of a specific 
element type is neglected in the analysis of the planar array 
since in normal usage it will have little effect for a planar 
array with a large number of elements. However, as with the 
linear array, the user of this software should be able to use 
information regarding the type of element and the orientation 
of the elements to narrow the range of the possible 6,,S and 
~,’S- 

The normalized array factor (AF) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [{Ref. 3: pp. 


260-263]: 


20 


. {| M . {| N 
sini — y,(6,) sins ¥ (8,9) 
AF,(9,) = eile ee) see | (4.1) 


mM sin| rE w sin| oe) 


where 
y,(0,0) = kd,sin(6)cos(@) +B, (radians) (4.2) 
p,(8,o) = kd,sin(8) sin(@) +B, (radians) (4.3) 
B, = -Kd,sin(8,) cos(@,) (radians) (4.4) 
B, = -kd,sin(6,) sin(9,) (radians) (4.5) 


Equation (4.1) is the normalized array factor that is a 
function of both 8 and @ where (4.2) and (4.3) are the array 
factor phase shifts in the x and y directions, respectively. 
(4.4) and (4.5) are the progressive phase shift between 
adjacent elements in the x and y directions, respectively, 


with the main beam directed along 0 = 9, and ¢ = @ In (4.2) 


°° 
through (4.5), k is the wavenumber (27/A). 

When the inter-element spacing between the elements (d, 
and d,) equals or is greater than the wavelength (A), multiple 
maxima of equal magnitude are formed. To avoid creating 
multiple maxima, the inter-element spacing must be less than 


a wavelength at the frequency of operation [Ref. 3: p. 263]. 


With the dependance of inter-element spacing on the 
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wavelength, planar arrays are characteristic narrowband 
antenna systems. 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r2>1.6A (meters) (4.6) 


SS (Md)? > (NGL)? (meters) (4.7) 


‘ 2M) se Cd.) 


7 (meters) (4.8) 


The length of an array is assumed to include a distance d/2 
beyond each end element [Ref. 5: p. 55]. Therefore, the 


maximum dimension of the planar array is: 


(MGUY Ss CNG). (meters) (4.9) 


The minimum distance to the far-field is found selecting the 
largest value of (4.6), (4.7), and (4.8). 
Directivity (D,) for planar arrays is determined 


approximately by [Ref. 7: p. 78]: 


AF AF; | 


dS 
M" 
p= 


a a a 4.10 
i i, “(aF,AF,)sin(@) dd db ( ) 


The maximum value of the numerator in (4.10) is unity at 0 = 
6 and @ = d, Since the array factor is normalized. The 


° 


denominator is the radiated power (P..,) of the antenna where 
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Piaq 1S Obtained by integrating over all angles around the 
antenna above the x-y plane (0 < ® < 7/2) [{Ref. 3: p. 28; Ref. 
7: p. 78). This assumes the array is always pointed toward 
one side or the other of the x-y plane. 

Effective isotropic radiated power (EIRP) for planar 
arrays is defined as the product of antenna gain and the total 
power radiated. EIRP is determined as [{Ref. 4: p. 62]: 


EIRP = PragD,  (W) (4.11) 


Assuming that the power delivered to planar arrays 
constructed of M x N isotropic radiating elements is 
equivalent to the radiated power, we estimate the radiation 
resistance (R,) of the array as [Ref. 3: p. 55]: 


= 2 ( Praa) 


oe FAL (Q) (4.12) 


Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a planar array (R,,) is 
assumed to be equivalent to the radiation resistance: 


(Q) (4.13) 


The gain (G) of a planar array cannot be calculated given 
the information assumed to be available, but an idealized 
value can be derived. Another factor affecting the gain is 
the assumption made when determining the total antenna 


efficiency (e€,). The total antenna efficiency for an antenna 


23 


is [Ref. 3: pp. 44-45): 


€, = €,€.7 = 1-(|- |? (dimensionless) (4.14) 
where 
bag = 1 (dimensionless) (4.15) 
e.=1-(|l |? (dimensionless) (4.16) 
Rin Zo ° . 
P= — = (dimensionless) (4.17) 
Rint Zo 


In (4.14), €,, is the conduction and dielectric efficiency, 
which is unity for an ideal lossless antenna. In (4.17), 2, 
is the characteristic impedance of the transmission line and 
r is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(4.15) when ideal lossless antenna is assumed. The ideal gain 
of a planar array in dimensionless and decibel quantities is 


then [Ref. 3: pp. 43-44]: 


G = €,D, (dimensionless) (4.18) 

G(dB) = 101log,, (€,D,) (dB) (4.19) 

The maximum effective aperture (A,,) for a planar array 

is derived by assuming polarization mismatches are negligible 


(this application is constructed using isotropic radiating 
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elements for which polarization loss factor (PLF) and antenna 


polarization are not calculated) [Ref. 3: p. 63]: 


Bo CANIS Bye (m2) (4.20) 
4u 


PLF =1 (dimensionless) (4.21) 


The effective height (h,,) of the planar array is now 


estimated by [{Ref. 6: p. 42]: 


hy .22 r**em (m) (4.22) 
No 
In (4.22), n, is the characteristic impedance of free space 


(1207). 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is in the x-z plane; 
the second is in the y-z plane; the third display generated is 
a plane perpendicular to the x-y plane when ¢ = 71/4; the last 
display is a plane perpendicular to the x-y plane when @¢ = @¢,. 
The last display is redundant when @ is either 0, 7/4, or 7/2. 

Results obtained with this application are the same as 
those found in [{Ref. 3: pp. 260-274]. Directivity was 
calculated using the approximate equation (4.10). An 
alternative approximation for directivity is given by either 


[Ref. 3: p. 272-273]: 
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D, = ~cos(@,) D,D (dimensionless) (4.23) 


° y 


or 


1? 


Ds = ————— (dimensionless) (4.24) 
Q, (rads?) 

In (4.23), D, and D, are the directivities of broadside linear 
arrays of length and number of elements Md,, M, and Nd,, N in 
the x and y directions, respectively. In (4.24), fy, is the 

beam solid angle. 
Table 4.1 is a comparison of directivity computed for a 
10x10 planar array. The reference directivity is derived 


using (4.24) while the calculated directivity is from (4.10). 


Table 4.1 Planar Array Data Comparison 


ANTENNA REFERENCE CALCULATED 


DIMENSIONS DIRECTIVITY (D,) DIRECTIVITY (D,) 


23.67 dB 24.07 dB 


-1/ (2/2) 


Figure 4.2 and Figure 4.3 are comparisons of array factor 


patterns for a 5x5 planar array with patterns from [Ref. 3] 


26 


and the calculated patterns. Figure 4.2 is a representation 
of the array factor patterns in the x-z and y-z planes, and 
the scale of this figure has been increased in order to 
increase the clarity of the pattern. Figure 4.3 is the array 
pattern in the plane at ¢@ = 45°, where the maximum beam is 
directed at ¢, = 45°. The array patterns produced by this 


application are identical to the patterns given in [Ref. 3]. 
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FIGURE 4.2 Array factor pattern of a 5x5 element, 
excited, planar array with d, = d A/2 m, 0, = 
45°, in the x-z and y-z planes 
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FIGURE 4.3 Array factor pattern of a 5x5 element, 
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in the plane @ = 45° 
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V. THE CIRCULAR ARRAY ANTENNA 


Circular arrays, or ring arrays, are antennas in which 
the radiating elements are placed on a circle with no elements 
positioned inside the circle. Circular antenna arrays with 
elements positioned inside the circle are planar arrays, and 
these antennas are addressed in Chapter IV. Circular array 
antenna are used for radio direction finding, air and space 
navigation, underground propagation, radar, and sonar systems 
{[Ref. 3: p. 274]. As with linear and planar arrays, circular 
arrays with uniformly spaced, equally excited elements are 
narrowband antenna systems due to the dependance of the number 
of elements (N) and the circumference of the circle to the 
wavelength (A) [Ref. 3: p. 277]. 


The following inputs are required to analyze a circular 


array: 
N = number of isotropic radiating elements 
a = radius of the circle 
f = frequency of interest “ 
6, = direction of main lobe at 6 = 6, 
@, = direction of main lobe at ¢ = ?, 


I. = antenna feed current * 
Z, = Characteristic feed impedance 
The first two inputs are data primarily obtained from visual 


data such as photographs, and the ~ indicates input parameters 
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that are either known or estimated. The planar array geometry 


is illustrated in Figure 5.1. 


FIGURE 5.1 Geometry of N-element circular array of 
isotropic point sources 
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This Mathcad application is written for N-element 
circular arrays with uniformly spaced, equally excited, 
isotropic radiating elements on the x-y plane. The 
polarization of the array is not calculated since isotropic 
point sources are used, and the antenna input resistance (R,,) 
is assumed to be equivalent to the radiation resistance (R,). 
This assumption results in the computation of ideal gain (G). 
As with the linear and planar arrays, the radiation pattern of 
a specific element type is neglected in the analysis of the 
circular array since in normal usage it will have little 
effect for a circular array with a large number of elements. 
However, as with both the linear and planar arrays, the user 
of this application should be able to use information 
regarding the type of element and the orientation of the 
elements to narrow the range of possible 6,’s and ¢,’s. 

The normalized array factor (AF) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [Ref. 3: pp. 


274-278; Ref. 4: 350-354): 


N 
AF(6,) = = e] {kasin®cos($-%,) + 4,1 (dimensionless) (5-1) 
n=1 
where 
®, = an()) n=1,2,...N (radians) (5.2) 
a, = ~kasin6,cos($,-®,) (radians) (5234 
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Equation (5.1) is the normalized array factor. Equation (5.2) 
is the angular position of the nth element on the x-y plane. 
Equation (5.3) is the phase excitation of the nth element with 
the position of the main beam directed at 8 = 6, and ¢ = 4¢,, 
and k is the wavenumber (27/A) and a is the radius of the 
array. 

An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r21.6A (meters) (5.4) 
Pe Sie) (meters) (5.5) 
re aiza)" (meters) (5.6) 


The minimum distance to the far-field is found by taking the 
maximum of (5.4), (5.4), and (5.6), where the diameter of the 
circle is taken to be the maximum dimension of the antenna. 
To determine the directivity (D,) for a circular array, 
the approximate equations of radiation intensity (U(6,¢)) for 


arrays and radiated power (P..,) are used. Radiation intensity 


rad 


is [Ref. 3: pp. 229-233]: 


U(6,o) = (AF)? (W/ solid angle) (5.7) 


The radiated power is [Ref. 3: pp. 28]: 


2n pn/2 : 
Pug] he U(8,) sin(6) dOdp  (w) (5.8) 
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In (5.8), P,,4 is computed by integrating over all angles above 
the x-y plane (0 < 8 s m/2) [Ref. 7: p. 78]. This assumes 
that the array has a pattern maxima above the x-y plane with 
no radiation below the x-y plane. The directivity for the 
circular array is [Ref. 3: pp. 229-233]: 


4nU. F ; 
ep max (dimensionless) (5.9) 


rad 
In (5.9), U,,, iS unity and occurs at 8 = 6, and @ = ¢, since 
the array factor is normalized. 

Effective isotropic radiated power (EIRP) is determined 


by [Ref. 4: p. 62]: 


EIRP = P,,,D, (W) (5.10) 


Assuming that the power delivered to a circular array of 
N isotropic, uniformly spaced radiating elements is equivalent 
to the radiated power, we can estimate then the radiation 


resistance (R,) of the array as [Ref. 3: p. 55]: 


Ro 2 ee (Q) (5.11) 


Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a circular array (R;,,) is 
estimated to be equivalent to the radiation resistance: 


(Q) (5.12) 
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Since the input resistance is estimated, the gain (G) of 
a circular array cannot be calculated, but an ideal value can 
be derived. As with linear and planar arrays, addressed in 
Chapter III and Chapter IV, respectively, the total antenna 


efficiency (€,) is (Ref. 3: pp. 44-45}: 


6 = £28.95 = (PP (dimensionless) (5.13) 
where 
Ses (dimensionless) (5.14) 
e. =1-(r/? (dimensionless) (5.15) 
Ro Bg ’ ’ 
r= —— (dimensionless) (5.16) 
Ry,+Z 


In (5.14), €.4 is the conduction and dielectric efficiency, 
which is unity for an ideal lossless antenna. In (5.16), Z, 
is the characteristic impedance of a transmission line, and [I 
is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(5.15) when an ideal lossless antenna is assumed, and the gain 
of a circular array is [Ref. 3: pp. 43-44}: 


G = €.D, (dimensionless) (5.17) 


G(aB) = 101og,, (€,D,) ( AaB) (5.18) 
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The maximum effective aperture (A,,) of a circular array 
is derived assuming polarization mismatches are negligible 
(this Mathcad application is constructed using isotropic 
radiating elements for which polarization loss factor (PLF) 


and antenna polarization are not calculated) [Ref. 3: p. 63]: 


2 6S 


2 
As. ri PLF- (m?) (5.19) 


PLF =1 (dimensionless) (5.20) 


The effective height (h,,) of a circular array is 


estimated by [Ref. 6: p. 42]: 


he, = 2 r‘tem (m) (5.21) 
No 
In (5.21), n, is the characteristic impedance of free space 


(1207). 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is the x-z plane; the 
second is the y~-z plane; the third is a plane perpendicular to 
the x-y plane when @¢ = 7/4; the last is a plane perpendicular 
to the x-y plane when @ = ¢,. The last display is redundant 
when @ is either 0, 7/4, or 1/2. 

The directivity of a uniform circular array approaches 
the number of elements (N) as the radius (a) of the array 


becomes very large [Ref. 3: p. 277]. Table 5.1 provides the 
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directivity of a uniform circular array as a function of the 


radius of the array. 


Table 5.1 Circular Array Directivity Data 


RADIUS (m) 


a 


Figures 5.2 and 5.3 provide comparisons of the array 


NUMBER OF DIRECTIVITY (D,) 


ELEMENTS 


——— 


10 13.70 dB 


radiation patterns computed for a 10-element, uniformly 
spaced, equally excited circular array with a radius of 1.59 
meters with those in [Ref. 3]. Figure 5.2 is the radiation 
pattern in the x-z plane, and the radiation pattern in the y-z 
plane is shown in Figure 5.3. For these patterns, 0, = ¢ = 
O°. The array radiation patterns generated by this 


application are identical to the patterns given in (Ref. 3]. 
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FIGURE S.2 Array factor pattern of a 10-element, uniformly 
spaced, equally excited circular array with a = 1.59 m, 


and 0, = ¢, = 0°, in the x-z plane 


38 


ear died 
o-f 

\ 

5) 

| 


o- 


“1-08 -0.6 -0.4 -0.2 


FIGURE 5.3 Array factor pattern of a 10-element, uniform 
spaced, equally excited circular array with a = 1.59 n, 
and 6, = ¢, = 0°, in the y-z plane 


to) 


39 


VI. THE FOLDED DIPOLE ANTENNA 


The folded dipole is a center-fed wire antenna that 
performs in a manner similar to that of a linear, center-fed 
thin-wire antenna with the added advantage of a more easily 
matched input impedance [Ref. 8: p. 7-36]. The folded dipole 
antenna consists of two parallel dipoles connected at the ends 
to form a narrow wire loop. The folded dipole antenna is 


illustrated in Figure 6.1. 


2a 


2b 


FIGURE 6.1 Geometry of folded dipole 
antenna 
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The folded dipole can be assumed to behave as a thin-wire 
dipole when the following are satisfied: the length of the 
antenna is much greater than the diameter of the feed 
conductor (L > 2a); the distance between the conductors is 
greater than the radius of the feed conductor (d > a); 
(2mad/A)? < 1; all conductors are assumed to be perfect [Ref. 
9: p. 172]; 

This Mathcad application calculates the parameters of a 
folded dipole by assuming the antenna behaves as a thin-wire 
dipole with an equivalent radius (a,). The current maximum is 
normalized to one amp. The equivalent radius is [Ref. 8: p. 


7-38]: 


2 2 
a, = exp a*in(a) + b*1n(b) + 2abln(a) (meters) (6.1) 


(a+b)? 
This application provides the characteristics of the antenna 
both in free space and positioned horizontally over the earth. 
Although the folded dipole positioned vertically over the 
earth occurs frequently in VHF/UHF communications, it is not 
addressed in this research. 


The following inputs are required to analysis a folded 


dipole: 
a = radius of feed conductor 
b = radius of second conductor 
ad = spacing between conductors 
L = length of antenna 


41 


h = height of antenna above ground plane 

f = frequency of interest “ 

€.g = conduction/dielectric efficiency of conductors “* 
€, = relative dielectric constant of ground plane “ 

o = conductivity of ground plane “* 


Yep = AGistance for far-field parameter calculations 
Oo, = incoming wave electric field unit vector for antenna 


in free space “ 


Oo, = unit polarization vector for antenna in free space” 

oh, = incoming wave electric field unit vector for 
horizontally positioned antenna “ 

6h, = antenna polarization direction for horizontally 


* 


positioned antenna 
The first five inputs are physical dimensions obtained through 
photographs, the * indicates input parameters that are either 


known or estimated. 


A. FOLDED DIPOLE IN FREE SPACE 

Assuming the folded dipole behaves as a thin-wire dipole 
with a, «< dad < L, we have the far-field electric field 
intensity of the antenna aligned with the z-axis as [Ref. 3: 
p. 120}: 


: kL kL 
Tea cos|— cos (8) me cos(=] 


HT Ie Bees Sin (6) 


In (6.2), n, is the intrinsic impedance of free space, equal 
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to 1207, and k = 27/A is the wavenumber. The normalized 
electric field is found by dividing (6.2) by the maximum value 
of the electric field which occurs at 6 = 7/2 for L less than 
about 1.25/. Clearly, the radiation pattern of the folded 
dipole is identical to that of a conventional dipole of length 
L. 

An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r21.6A (meters) (6.3) 
r>s5sL (meters) (6.4) 
C2 ai (meters) (6.5) 


The minimum distance to the far-field is found by comparing 
the values of (6.3), (6.4), and (6.5) and selecting the 
maximum value. The distance between the conductors is not 
taken into account in determining the minimum distance to the 
far-field since it is assumed that the folded dipole has an 
over all effective radius (a,). 

Directivity (D,) for the folded dipole is found by 
determining the radiation intensity (U(8)) and radiated power 
(P.aa)- Since the radiation pattern is symmetric in the x-y 
plane, U(@) and P,,, are functions of @ only. Radiation 


intensity is [Ref. 3: p. 28]: 
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2 
U(6) = SH |B(8) | (W/ solid angle) (6.6) 


The radiated power is [Ref. 3: pp. 28]: 


Prag = 2n{"U(8) sin(@) cB (W) (6.7) 
ie) 


Directivity for a folded dipole is [{Ref. 3: pp. 29]: 


4nU : 
5 = (dimensionless) (6.8) 


rad 
where U,,, = U(1/2). 

Effective isotropic radiated power (EIRP) is determined 
by {Ref. 4: p. 62]: 


EIRP = P.agD,  (W) (6.9) 


The folded dipole is unique in its capability to act as 
a step up impedance transformer when the antenna is resonant. 
The antenna effectively operates as an unbalanced transmission 
line, and it can be analyzed by decomposing the current into 
two distinct modes: a transmission line mode [Figure 6.2] and 
an antenna mode [Figure 6.3]. {Ref. 3: p. 341; Ref. 4: pp. 


205-207] 
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FIGURE 6.2 Transmission FIGURE 6.3 Antenna mode 
line mode 


The input impedance (Z;,,) at the feed for the folded 
dipole is given in (6.10), and the radiation resistance (R,,) 
for the antenna is the real component of (2Z,,)) [Ref. 3: 
p.342]: 


ee ec a (6.10) 
Be ae Oe 
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In (6.10), Z, is the input impedance of the transmission line 


mode and is obtained from the impedance transfer equation 


= jZ tan} (Q) (6.11) 


where Z, in (6.11) is the characteristic impedance of a two- 
wire transmission line and is given in (6.12) [{Ref. 3: pp. 


341-342; Ref. 8: p. 7-39]: 


Zo = 7? cosh” {45 (Q) (6.12) 


2fab 
The input impedance of a linear dipole of length L and 
equivalent radius a, is calculated by [Ref. 3: pp. 124, 127, 


294, and 342]: 


Za = Ry + JjXq (Q) (6.13) 
R 
Ry 2 (Q) 
sin <2) (6.14) 
X 
fee, SC) 6.15 
sin(<*) 


R, = ely + In(kL) - ¢,(KL)] 


+ Belg sin (kL) (8, (2KL) - 28,(kL)]| (6.16) 
ois 2 


Sealy cos (KL) (r + In| =) + C,(2kL) -20;(kL)} 
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3 


X, = 72 [25;( kL) +cos (kL) [25;(kL) -5;(2kL)]] 
(6.17) 
2 
sale sin( kt) [2¢,(kE) -¢,(2kL) ~cf 22 )) 
4Tm L 


In (6.16) and (6.17), C,(x) and S,(x) are the cosine and sine 


integrals, respectively, [Ref. 3: pp. 124 and 743-746]: 


C, (x) = [7 S28) ae = yein(x) -y (6.18) 
© n=1 : 

: x sin(t) 2 ie ( -1) 271 x (20-1) 6.19 

a I, =o 2 (2n-1)(2n-1)! 


and y 0.57721 (Euler’s constant). 

Mathcad computes summations more quickly than integrals. 
In this application, the summation forms of (6.18) and (6.19) 
are used with the index n ranging from 1 to 50 to ensure 
proper convergence. The maximum value of the argument of 
either C,(x) or S,(x) is x = 4mL/d, and for L = A the 
succeeding terms in the summation rapidly decrease for n > 20. 

The gain (G) of a folded dipole is obtained as the 
product of the antenna reflection efficiency (€,) and the 
directivity of (6.8) where the antenna reflection efficiency 


is [Ref. 3: pp. 43-45]: 


€, = €,€04 = Egg(1-|F|?) (dimensionless) (6.20) 
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where 


e = 1-(|P |? (dimensionless) (6.21) 
Rind ; ; 

r= —“—_ (dimensionless) (6.22) 
Rin + Zo 


In (6.21) e. is the mismatch efficiency and e is the 


r cd 


conduction and dielectric efficiency in (6.22). In (6.22), 2, 


is the characteristic impedance of the transmission line given 


in (6.12). Therefore, the gain of the antenna is: 
G = ¢€,D, (dimensionless) (6.23) 
G(dB) = 101log,, (€,D,) (dB) (6.24) 


The folded dipole in free space aligned on the z-axis is 
horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave’s unit 
polarization vector (9,), the polarization loss factor (PLF) 


is [Ref. 3: pp. 48-53]: 


PLF = |G,- 63 |? (dimensionless) (6.25) 


The maximum effective aperture (A,,) of a folded dipole 


is (Ref. 3: p.63]: 


A - SAY pir (m? ) (6.26) 
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The effective height (h,,) of the folded dipole is (Ref. 


6: p. 42]: 


h.=2 r**em (m) (6.27) 


The bandwidth of the folded dipole antenna is dependant 
on the length of the wires and has better bandwidth 
characteristics than a conventional dipole of the same length. 
It can be assumed that the bandwidth of the folded dipole is 
essentially the same as that of a conventional dipole of an 
radius a, (a < a, < d/2) [{Ref. 3: p. 346]. The bandwidth is 
increased when the distance between the conductors (d) is 
increased and/or when the conductor radii a and b are 
increased (Ref. 10: p. 190]. 

Tables 6.1 and 6.2 are comparisons of measured and 
calculated input impedance data for the folded dipole as a 
function of its electrical length. The radii of the 
conductors (a, b) are the same (0.0005 m) while the spacing 
between the conductors (d) is held constant at 0.00625 m and 


4 = 1m (Ref. 8: pp. 7-40 - 7-41}. 
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TABLE 6.1 Folded Dipole Antenna Data Comparison 


CALCULATED 


ANTENNA MEASURED 


LENGTH (L) 


Re(Z,,) 2 Re(Z,,) 2 


xx 


A/a 6533 


2592 
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TABLE 6.2 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED CALCULATED 


LENGTH (L) Im(Zi,) 2 Im(Z,,) 2 


To further test the validity of the Mathcad application, 
a numerical antenna analysis program was also used. The 
program used was ELNEC, a powerful, easy-to-use program for 
modeling and analyzing a wide variety of antenna types 
including ground effects and parasitic structures. The 
fundamental computation portion of ELNEC is the same as 
MININEC (Version 3), that was developed by the Naval Ocean 
Systems Center [{Ref. 11: pp. 5-7]. 

Table 6.3 is a comparison of the gain for a folded dipole 
in free space with different dimensions [Ref. 8: p. 3-29]. 


For each case, the radii of the conductors (0.0005 m) and the 
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distance between the conductors (12.5a m) are the same. As 
can be seen, the comparison between the gains obtained from 


various sources is very good. 


TABLE 6.3 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED ELNEC Mathcad 


DIMENSION Gain (G) Gain (G) Gain (G) 


1.70 GB 


The computed parameters in Tables 6.1, 6.2, and 6.3 are 
essentially the same as those measured or computed with ELNEC 
with exception of the ones marked with "". These discrepancies 
result in using the transmission line model and equivalent 
radius to analysis the folded dipole in this application. 

As previously mentioned, the folded dipole is a center- 
fed wire antenna that performs in a manner similar to a 
linear, center-fed, thin-wire antenna. Therefore, the 
radiation pattern of a folded dipole is also similar to that 
of a thin-wire antenna. The folded dipole radiation pattern 


in free space is illustrated in Figure 6.4. 
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FIGURE 6.4 Radiation Pattern folded dipole in free space 
with L = A/2 m, a = b = 0.0005 m, and a, = 0.004 m 


Figure 6.5 is a comparison of the radiation pattern of a 
folded dipole in free space obtained using ELNEC and the 
Mathcad application. The patterns are not identical, but are 
Similar and oriented in the same direction as expected. The 
difference can be attributed to the fact that ELNEC takes into 


account mutual coupling between the conductors while the 
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Mathcad application models the folded dipole as a thin-wire, 


center-fed dipole. 
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Figure 6.5 Comparison of radiation patterns using ELNEC 
and Mathcad of a Folded dipole in free space with L = A/2 
m, a = b = 0.0005 m, and a, = 0.004 m 
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B. FOLDED DIPOLE POSITIONED HORIZONTALLY OVER THE EARTH 
Through the use of image theory, the far-field electric 
field intensity of a folded dipole oriented horizontally above 
a flat earth can be obtained. The @ component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 


Space is (Ref. 3: pp. 143-144]: 


Klos ; kL 
-jkr¢| COS{/ —— sin (6) sin(p)| -cos|— 
£0 (8, ) spp eee aoe ee) (v/m) (8-28) 
2h Ere V1 -sin? (6) sin) 


From image theory, the 8 component of the electric field 


over ground is [Ref. 4: pp. 229-235]: 


E0,, (8) =E0(0,) (er icos (©) a Per eaeee (V/m) (6.29) 


In (6.28) and (6.29), E,(0,¢) is valid only above the ground 
(O < 80 < mw/2) (Ref. 3: pp. 135-142}. In (6.29), h is the 
distance from the earth to the antenna and the vertical 


reflection coefficient (I) of the ground is (Ref. 4: pp. 231- 


233]: 
/ - fe’,- sin? (6) 
| Ras = £' cos (8) - ye',- sint (8) (dimensionless) (6.30) 
e’ cos(6) + Je’,- sin? (6) 
where 
boos _ o , , 
es bpe 7 Dn fe, (dimensionless) (6.31) 
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In (6.30) and (6.31), (€’,) is the relative complex effective 
dielectric constant of the ground, (€,) is the relative 
dielectric constant of the ground, and (ga) is the conductivity 
of the ground. For a perfect ground plane, 0 > ~ and Yl, = 1. 

The # component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 
135-142}: 


-jkr ¢¢ cos| Esin (6) sin(>)} -cos(— 


| (V/m) 


I 
Eb(0,) =jn,—2 


Or. 1 -sin? (6) sin*(6) 


(6.32) 


Now from image theory, the ® component of the electric field 


over earth is [Ref. 4: pp. 229-235]: 


Eb, (8) =£ (0, ) (eF#heo8(®) 4 Pie-Fkheos(®) = (v/m) (6.33) 


In (6.33), (T,) is the horizontal reflection coefficient of 


the ground and is [Ref. 4: pp. 231-232]: 


— _ } 2 
Gee ea ee eieeenieas: 06.50) 


cos (6) + /e7,- sin? (6) 


When o + ©, T, = -1l. 

The radiation intensity (U,(8,¢)) of a folded dipole over 
a earth must take into account both the 6 and # components of 
the electric field. Therefore, the radiation intensity is 


[Ref. 3: p. 28]: 
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pare 
3! ap () |? + |Bb,,(8) [7] (W/ solid angle) (6.35) 


U,, (8) = 


The radiated power (Ph,,,) follows (6.7) except that the field 


is present only above the earth. Thus, 


Phyag = 20 f 7 U,(8) sin(8) d8 (Hi) (6.36) 


Directivity (Dh,) and effective isotropic radiated power 
(EIRPh) for a horizontal folded dipole over the earth are 
identical to (6.8) and (6.9), respectively, with the results 
of (6.36) inserted. 

Input resistance (Rh,,) for the folded dipole over the 
earth plane is found by (where I, in (6.32) is normalized to 


one amp) [Ref. 3: p. 124]: 


in sin £2) (6.37) 


where the radiation resistance (Rh,) is, 


2 Ph 
Rh, = ——22 = (Q) (6.38) 

2 

IZ3| 
The voltage reflection coefficient (ITH), gain (Gh), 
maximum effective aperture (Ah,,), and maximum effective height 
(hh,,) for a horizontal folded dipole over the earth are 


determined using the equations for those parameters for a 


folded dipole in free space, as listed above. 
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The Polarization Loss Factor (PLFh) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (8h,) is estimated. The antenna 
polarization vector (oh,) is then calculated, and estimating 
the unit polarization vector of the incoming wave (oh,), the 
PLF is found from (6.25). 

Table 6.4 is a comparison of measured [Ref. 8: p. 3-29], 
ELNEC, and Mathcad software gains for a horizontal folded 
dipole over a reflecting sheet. The radii of the conductors 
(0.0067 m) are the same. As can be seen, the Mathcad software 


gain is similar to the measured and ELNEC gains. 


TABLE 6.4 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED ELNEC Mathcad 


DIMENSION Gain (G) Gain (G) Gain (G) 


A/2 m 


d = L/25.5 m 


A/8 m 
eee ee 


No input impedance data has been found to compare a 
folded dipole oriented horizontally over the earth with that 
computed by this Mathcad application. 

Figures 6.6 and 6.7 are the H-plane and E-plane radiation 


patterns of a folded dipole oriented horizontally over the 
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earth, respectively. Since the folded dipole is similar to a 
thin-wire center-fed dipole of length L, the illustrations 


below are similar to those of a dipole of length L. 


FIGURE 6.6 H-plane radiation pattern of a folded dipole 
with L = A/2 m, a = b = 0.0005 m, ad = 0.00625 m, h = 0.5 
m, 0 = 0.01 S/m, e, = 15 


r 
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FIGURE 6.7 E-plane radiation pattern of a folded dipole 
with L = A/2 m, a = b = 0.0005 m, da = 0.00625 m, h = 0.5 
m, 0 = 0.01 S/m, ¢, = 15 


Figures 6.8 and 6.9 are comparisons of the H-plane and E- 
plane radiation patterns between those obtained with ELNEC and 
Mathcad software of a folded dipole oriented horizontally over 
the earth, respectively. As expected, the patterns are 
Similar. As previosly mentioned, the difference can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the conductors while the Mathcad application 


is modeled after a thin-wire, center-fed dipole. 
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Figure 6.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a folded dipole with L = A/2 m, a 

b = 0.0005 m, dad = 0.00625 m, h = 0.5 m, o = 0.01 S/m, € 
15 


r 


0.4 
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Figure 6.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a folded dipole with L = A/2 m, a = 
b = 0.0005 m, a = 0.00625 m, h = 0.5 m, o = 0.01 S/m, €, = 
15 
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VII. THE CAGED DIPOLE ANTENNA 


The caged dipole antenna is used in communications for 
its broadband antenna characteristics. Its performance is 
similar to that of a thick cylindrical dipole. Since thick 
cylindrical dipoles are usually too clumsy to use except at 
very short wavelengths, thick dipoles are replaced with a 
caged dipole which consists of a number of thin wire dipoles 
arranged to form a circular cross section. The antenna is 
usually supported by a central conductor and support members 
that provide rigidity but contribute little to the electric 
properties of the caged dipole. The caged dipole (without the 
central conductor and the support members) is illustrated in 


Figure 7.1 (Ref. 12: pp. 172-174]. 


FIGURE 7.1 Geometry of caged dipole antenna 
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This Mathcad application models the caged dipole antenna 
as a circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles. This approach is taken to account 
for the number of conductors that are present in the caged 
dipole. As with the arrays previously considered, mutual 
coupling between the dipoles is neglected. In this manner, 
this application simulates the characteristics of a thick 
cylindrical dipole. The thin-wire dipoles are modeled after 
finite length dipoles in free space that have their length 
much greater than their diameter (L > wire diameter), and the 
dipoles are assumed to be perfect conductors. Since the 
central conductor and the support wires have little effect on 
the electrical properties of a caged dipole, mutual coupling 
from these components are not taken into account in this 
application. As seen in Figure 7.1, the inner ends of a caged 
dipole, the feed points, maybe coned in order to reduce base 
capacity [Ref. 12: pp. 172-174}. This effect is also not 
addressed. 

This application provides the characteristics of the 
antenna in free space as well as and oriented both vertically 


and horizontally over the earth. 


A. CAGED DIPOLE IN FREE SPACE 
The following inputs are required to analyze a caged 


dipole in free space: 


63 


N = number of conductors 
a = radius of antenna 


L = length of antenna 


Au 


f = frequency of interest 

I, = antenna feed current * 

Z, = input impedance * 

£.4 = conduction/dielectric efficiency of conductors “ 
Yr, = Adistance for far~field parameter calculations * 
Oo, = incoming wave electric field unit vector for 


Ld 


antenna 


ome unit polarization vector for antenna * 
The first three inputs are physical dimensions obtained 


ve 


primarily from photographs, and the indicates input 
parameters that are either known or estimated. 

Assuming the caged dipole behaves as a circular array 
with equally excited (I,), uniformly spaced elements in the x- 


y plane with a radius of (a), we have the array factor 


(AF(6,¢)) [Ref. 3: pp. 274-278]: 


N 
AF (6,6) = > yO cl ae A (dimensionless) (7.1) 


n=1 


where the angular position in the nth element in the x-y plane 


is: 


®, = an(a)en SOs, ck Se, (radians) (7.2) 
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Equation (7.1) is derived from: 


N 
AF (0,6) = SY T,e7 {X28!n0008(%- Pn) +e) (dimensionless) (7-3) 
n=1 
where the phase excitation of the nth element is: 


a, = 0 = -kKasin6,cos(o,-®,) (radians) (7.4) 


In (7.3), and the phase excitation of the nth element (a,) is 
zero as the result of desiring the main beam of the array to 
be directed perpendicular to the length of the antenna (8, = 
m/2). In (7.1), k is the wavenumber (27/A). 

Assuming each conductor of a caged dipole behaves as a 
thin-wire dipole, we have the far-field electric field 
intensity of a single conductor aligned with the z-axis as 


{[Ref. 3: p. 120): 


Te dk cos( <= cos (8)) - cos{ =} (7.5) 
E(8) =jn,— (V/m) : 
2UT ep sin (6) 
In (7.5), n, = 1207 and is the intrinsic impedance of free 


space. 

As previously discussed in Chapter III, the total 
electric field of an array constructed with identical elements 
is the product of the field of a single element and the array 
factor as a consequence of the principal of pattern 


multiplication (Ref. 3: p. 207). Therefore, the total 
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electric field of a caged dipole in free space is approximated 
as: 


E.(6,@) = £(8) AF(8,0) ( V/m) (7.6) 


An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r2>1.6Aa (meters) (7.7) 


r2 5VL? + (2a)? (meters) (7.8) 


2(L2 + (2a)?) 


FE .2 ry 


(meters) (7.9) 


The minimum distance to the far-field is found by taking the 
maximum of (7.7) through (7.9). In (7.8) and (7.9), the 
maXimum dimension of the antenna is the diagonal length of the 
array. 

Directivity (D,) for a caged dipole is found by 
determining the radiation intensity (U(8,¢)) and radiated 


power (P,,4)- Radiation intensity is [Ref. 3: p. 28]: 


2 
U(8,o) = a |E.(8,@) |? (W/ solid angle) (7.10) 


The radiated power is [Ref. 3: pp. 28]: 


2n rn , 
Pag. iB U(8,@) sin(6) dddp (Ww) (7.11) 
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Directivity for a caged dipole is [Ref. 3: pp. 29]: 


4nU , ; 
ee a (dimensionless) (7.12) 


rad 


where U,,, = U(7/2,0). 
Effective isotropic radiated power (EIRP) is determined 


by [Ref. 4: p. 62): 


EIRP = P,.gD,  (W) (7.13) 


Input resistance (R,,) for a caged dipole in free space 


is found from (Ref. 3: p. 124}: 


R 
Rk Se a es Q 
in sini =| ( ) (7.14) 
2 
where the radiation resistance (R,) is: 
AF 
R, = zane (Q) (7.15) 
|Z,,|? 


The gain (G) of a caged dipole is obtained as the product 
of the antenna reflection efficiency (e€,) and the directivity, 


where the antenna reflection efficiency is ({Ref. 3: pp. 43- 


45): 


€, = €,€o9 = Ecg(l-|F|?) (dimensionless) (7.16) 


where 


e- = 1+ |T |? (dimensionless) (7.17) 
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T= ——< (dimensionless) (7.18) 


The parameter e, is the mismatch efficiency, e,,; is the 
conduction and dielectric efficiency which is unity for a 
lossless conductor, and Z, is the characteristic impedance of 
a transmission line. Therefore, the gain of the antenna is: 


G=et,D, (dimensionless) (7.19) 
G(aB) = 101log,, (€,D,) (dB) (7.20) 


A caged dipole in free space aligned parallel to the z- 
axis is horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave’s unit 


polarization vector (o,), the polarization loss factor (PLF) 
is [Ref. 3: pp. 48-53]: 


PLF = |G,° G;|? (dimensionless) (7.21) 


The maximum effective aperture (A,,) of a caged dipole is 


(Ref. 3: p.63]: 


A. = SUA)" PLP (m?) (7.22) 


The effective height (h,,) of the caged dipole is [Ref. 


6: p. 42]: 


h..=2 r/*em (m) (7.23) 


The bandwidth of a thick cylindrical dipole antenna is 
dependent on the _ ratio of the length of the antenna to the 
diameter of the antenna (L/2a). When L/2a *% 5000 the 
acceptable bandwidth of a thick cylindrical dipole is about 3% 
of the center frequency, and when L/2a * 260 the bandwidth is 
about 30%. Hence, if we assume a linear relationship between 
bandwidth and L/2a, and since the caged dipole’s performance 
is similar to that of a thick cylindrical dipole, then the 
bandwidth of a caged dipole is [Ref. 3: p. 333): 


BW = fyic, = f (Hz) (7.24) 


low 


where the upper and lower frequencies are: 


Enign = fq + F( 35) f. (Hz) (7.25) 
a Psa} f (Hz) (7.26) 
low Cc 2a c 


In (7.25) and (7.25), f, is the center frequency and F(L/2a) 


c 
is the linear relationship function between antenna length and 
diameter. This relationship is obtained by solving y = mx + 
f, for m using x, = 5000, y, = 3% and x, = 260, y, = 30%. The 


relationship becomes: 


F[35]=-5-696-10° (35} 10-3248 (dimensionless) (7.27) 
a 
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Table 7.1 is a comparison of the gain of a thick 
cylindrical dipole with that of a caged dipole (N = 8) in free 
space computed both by ELNEC and Mathcad software [Ref. 8: p. 
2-38]. As previously mentioned in Chapter VII, ELNEC is an 
easy-to-use numerical analysis program for modeling and 
analyzing different antenna types [Ref. 11: pp. 5-7]. As can 
be seen, the comparison between the measured, ELNEC, and 


Mathcad application values are similar. 


TABLE 7.1 Caged Dipole Antenna Data Comparison 


ANTENNA MEASURED ELNEC Mathcad 


DIMENSIONS Gain (G) Gain (G) Gain (G) 


Figure 7.2 is the radiation pattern of a caged dipole (N 
= 8) in free space. It is essentially equivalent to that of 
a center-fed, thick cylindrical dipole having the same 


conductor length and diameter [Ref. 8: p. 2-38]: 
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FIGURE 7.2 Radiation pattern of a caged dipole in free 
space with L = A/2 m, N = 8, and a = 0.125 m 


Figure 7.3 is a comparison of the radiation pattern of a 
caged dipole (N = 8) obtained using ELNEC and Mathcad. The 
patterns are not identical, but are similar and oriented in 
the same direction as expected. The difference can be 


attributed to the fact that ELNEC takes into account mutual 
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coupling between elements while the Mathcad application 


ignores these effects. 
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FIGURE 7.3 Comparison of radiation patterns using ELNEC 
and Mathcad of a caged dipole in free space with L = A/2 
m, N = 8, and a = 0.125 m 
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B. CAGED DIPOLE ORIENTED VERTICALLY OVER EARTH 

The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented vertically 
over the earth. The height of the antenna above ground is a 


physical dimension obtained through photographs: 


h = antenna height above ground 

€, = relative dielectric constant of earth 

o = conductivity of earth 

ov, = incoming wave electric field unit vector for 


vertical antenna 

6v, = antenna polarization direction for vertically 

orientated antenna 

As with the caged dipole in free space, the analysis of 
a caged dipole positioned vertically above the ground does not 
take into account the central conductor and support members 
Since they produce negligible effects on the antenna’s 
electrical properties. 

From image theory and pattern multiplication for arrays, 
the total far-field electric field intensity of a caged dipole 
oriented vertically above a flat earth is (Ref. 4: pp. 229- 
235): 


Ev, (6,6) =Ev(8) AFV(6, o) (e7kRo0s'®) + Le Ferer (V/m) (7.28) 


In (7.28), Ev(6) is the far-field electric field intensity of 
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a finite length dipole in free space which is given by (7.5), 
and AFv(8,¢) is the array factor for a circular array which is 
(7.1). The total electric field intensity is valid only above 
the ground (0 < @ < 7/2) [Ref. 3: pp. 135-142]. The distance 
from the earth to the center of the antenna is (h), and the 
vertical reflection coefficient (I,) of the ground is [Ref. 4: 
pp. 231- 233]: 


e’ cos (8) - /e’,-sin? (6) 


e’ cos (6) + fe’, -sin? (6) 


r,(6) = (dimensionless) (7.29) 


where 


‘ 0 


e/ = od eae (dimensionless) (7.30) 
oO 


I I 


In (7.29) and (7.30), (€’,) is the relative complex effective 
dielectric constant of the ground, (¢€,) is the relative 
dielectric constant of the ground, (co) is the conductivity of 
the ground, and ¢, is the permittivity of free space (1/(367) 
x 10°). For a perfect ground plane, o + ~ and TI, = 1. 

The radiation intensity (U,(8,¢)) of a caged dipole over 


earth is [Ref. 3: p. 28]: 


2 


u,(8,o) = —££ 
21, 


[lZv, (8,9) |?] (W/ solid angle) (7.31) 


The radiated power (Pv,,3) is present only above the earth. 


Thus, [Ref. 3: p. 28]: 
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Pvraq = { [2 Uv(8,o) sin(®) ddd = (w) (7.32) 
0 ie} 


Directivity (Dv,) and effective isotropic radiated power 
(EIRPv) for a vertically oriented caged dipole over the earth 
are given by (7.12) and (7.13), respectively, with the results 
of (7.32) inserted. UVmax in (7.12) is simply the maximum 
value of (7.31) which, in Mathcad, is found by constructing 
(7.31) as a vector and then finding the maximum value of the 
vector. 

Input resistance (Rv,,), voltage reflection coefficient 
(TV), gain (Gh), maximum effective aperture (Av,,), and maximum 
effective height (hv,,) for a vertical caged dipole over the 
earth are determined using the equations for those parameters 
for a caged dipole in free space. 

The Polarization Loss Factor (PLFv) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (6v,) is estimated. The antenna 
polarization vector (ov,) is then calculated and, estimating 
the unit polarization vector of the incoming wave (ov,), PLF 
is found from (7.21). 

Table 7.2 is a comparison of the gain of a vertical caged 
dipole over a flat earth obtained with both ELNEC and the 
Mathcad application. As can be seen, the gains are Similar. 
As previously mentioned, the difference in the gains can be 


attributed to the fact that ELNEC takes into account mutual 
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coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 


and approximates ground effects. 


TABLE 7.2 Caged Dipole Antenna Gain Comparison 


ANTENNA ELNEC Mathcad 


DIMENSIONS Gain (G) Gain (G) 


Since the caged dipole behaves in a manner similar to 
that of a thin dipole except for bandwidth characteristics, 
the radiation pattern of a vertically oriented caged dipole 
over ground will resemble that of a thin vertical dipole over 
ground. Figure 7.4 is a plot of the radiation pattern of a 
caged dipole over earth. There is a very close resemblance to 
the radiation pattern of a thin vertical dipole over ground 
[Ref. 4: p. 234]. Figure 7.5 is a comparison of radiation 
patterns computed by ELNEC and Mathcad software. As expected, 


the patterns are similar. 
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FIGURE 7.4 E-plane radiation pattern of a caged dipole 
with L = 4/2 m, N= 8, a = 0.125 m, h = 0.25 m, o = 0.01 
S/m, and e, = 15 


-— Mothcad Pattern, ... ELNEC Pottern 


FIGURE 7.5 E-plane radiation pattern comparison between 
ELNEC and Mathcad of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m, h = 0.25 m, o = 0.01 S/m, and ¢€, = 15 
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C. CAGED DIPOLE ORIENTED HORIZONTALLY OVER EARTH 

The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented horizontally 
over the earth: 

h = antenna height above ground 

€, = relative dielectric constant of earth 

o = conductivity of earth 

oh, = incoming wave electric field unit vector for 

vertical antenna 

6h, = antenna polarization direction for vertically 

orientated antenna 

As with the evaluations of a caged dipole both in free 
space and oriented vertically over the earth, the analysis of 
a caged dipole positioned horizontally above the ground does 
not take into account the central conductor and support 
members since they produce negligible effects on the antenna’s 
electrical properties. 

Through the use of image theory, the far-field electric 
field intensity of a caged dipole oriented horizontally above 
a flat earth can be obtained. The 8 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 


Space is [Ref. 3: pp. 143-144]: 
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0 EE 


es gor 


come P 
-jkr¢¢| COS( —— Sin (8) sin(o)| -cosj —— 
£0h (8,6) =jn,—2~ G28 ge EO) 0a 


2 V1 -sin? (6) sin%() 
The 8 component of the total electric field over ground 


is [Ref. 4: pp. 229-235]: 


EOh, (8, ) = £0h(8,) AFA(O, ) (eF*h-05(9) _ Te -Fkheos(®)) — (y/m) (7.34) 


The array factor (AFh(6,q¢)) in (7.34) for a horizontal array 


is determined by [Ref. 3: pp. 274 and 776]: 


N 
AFhH(8,o) = D> 1, e7%*8*78S*neSin%, (dimensionless) (7-35) 
n=1 


In (7.33) through (7.35), the electric field (E6h(9,¢)), the 
total electric field (E6h,(8,¢)), and array factor (AFh(6,¢)) 
are valid only above the ground (0 < 8 < 7/2) (Ref. 3: pp. 
135-142]. In (7.34), (h) is the distance from the earth to 
the center of the antenna, and the vertical reflection 
coefficient (T,) of the ground is given by (7.29). 

The © component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 
135-142]: 


pe Jkt cos{ “sin (8) sin(b)}-cos(=5 
Eoh(6,o) =jn,— —— 


OE ge 1-sin? (6) sin*(6) 


From image theory, the ® component of the total electric field 


over earth is [Ref. 4: pp. 229-235]: 
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Ebh, (0,6) =Bbh(8,) AFH(8,o) (eF*ors) + Tie Jkneos())  (v/m) (7.37) 


In (7.37), (AFh(8,¢)) is given by (7.35), and (T,) is the 
horizontal reflection coefficient of the ground and is [Ref. 
4: pp. 231-232]: 


= / 7 es se 
I, COSCO) rye ee AO) (dimensionless) (7.38) 
cos (@) + fe’,- sin? (6) 


When o - ©, Ty, = -1. 

The radiation intensity (U,(8,¢)) of a caged dipole over 
a earth must take into account both the 8 and ¢# components of 
the electric field. Therefore, the radiation intensity is 


[Ref. 3: p. 28]: 


2 
Uh (8,6) = [|20h, (0,6) |? +|Bbh, (0,6) |?] (W/ solid angle) (7.39) 


The radiated power (Ph,,,;) is present only above the earth. 


Thus, 


2n pn/2 : 
PRrag = f i Uh(0,@) sin(®) dOdb (wv) (7.40) 


The remaining parameters (i.e., Dh,, EIRPh, Rh,,, etc.) 
for a horizontal caged dipole over the earth are calculated 
using the equations for these parameters listed previously. 


Table 7.3 is a comparison of the gain of a horizontal 


caged dipole over a flat earth obtained with ELNEC and the 
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Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in gains can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 


and approximates ground effects. 


TABLE 7.3 Caged Dipole Antenna Gain Comparison 


ANTENNA ELNEC Mathcad 


DIMENSIONS Gain (G) Gain (G) 


0.01 S/m 


= 15 


r 


As previously mentioned with regard to the vertically 
oriented caged dipole over the earth, the caged dipole behaves 
in a manner similar to that of a thin cylindrical dipole. 
Therefore, the radiation patterns of a horizontally oriented 
caged dipole are similar to those of horizontal thin 
cylindrical dipole which in turn resembles that of a thin 


dipole. Figures 7.6 and 7.7 illustrate the H-plane and E- 
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plane patterns, respectively, of a caged dipole over the earth 
(Ref. 4: pp. 231-237). The radiation patterns of the 
horizontally orientated caged dipole over a flat earth are 
very similar to those of a horizontal thin dipole over earth. 
Figures 7.8 and 7.9 compare the radiation patterns between 
those computed with ELNEC and those computed by Mathcad 


software. As expected, the patterns are similar. 


FIGURE 7.6 H-plane radiation pattern of a caged dipole 
with L = A/2 m, N = 8, a = 0.125 m, h = 1.0m, o = 0.01 
S/m, and e€, = 15 
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FIGURE 7.7 E-plane radiation pattern of a caged dipole 
with L = A/2 m, N = 8, a = 0.125 m, h = 1.0 m, o = 0.01 
S/m, and e, = 15 


0.4 


Pattern 


FIGURE 7.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m, h = 1.0 m, o = 0.01 S/m, and e, = 15 
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FIGURE 7.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m, h= 1.0m, o = 0.01 S/m, and ¢, = 15 
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VIII. PARABOLIC REFLECTORS 


Many communications and radar operations require very 
large values of antenna gain (G) that are difficult to achieve 
with a single device. One of the most popular ways to 
construct a very high gain antenna is to reflect the 
electromagnetic energy of a relatively small antenna off a 
metallic parabolic dish in the far-field of the feed device 
(Ref. 13: p. 12-2]. If the source antenna is located at the 
focus of the parabola and if the radius (a) of the parabolic 
dish periphery is large with respect to wavelength (A), then 
the reflector may be approximated as a uniform phase aperture 
antenna. Gain in excess of 30 dB is common for this type of 
structure [Ref. 4: pp. 423-424]. The geometry of a parabolic 
reflector is illustrated in Figure 8.1, where the parameters 
include the focal distance (f), the distance (p) from the 
origin to any point (P) on the reflector surface, and the 
distance (r’) from the origin to the projection of point (P) 
onto the z = O plane. The primed angles in Figure 8.1 
correspond to the parameters associated with the feed horn. 
The Mathcad application user should note that the focal 
distance should exceed the minimum distance to far-field for 


the feed antenna. 
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FEED HORN 


REFLECTOR 


FIGURE 8.1 Parabolic Reflector Geometry 


The Mathcad parabolic reflector applications are unique 
in that the software is not written as a stand alone package. 
Rather, the reflector software is run as an addendum to the 
programs that already exist for various antennas [Ref. 2: pp. 
4-15, 67-114}. Four types of feeds are analyzed: helical, 


planar equiangular spiral, pyramidal horn, and conical horn. 
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The following additional inputs are required to analyze 
parabolic reflectors: 
foc = reflector focal length 


arep = radius of the mouth of the reflector 


E, = electric field scale factor “ 

€, = blockage efficiency “ 

€., = spar efficiency * 

ero = Secondary far field observation distance “ 
X,Y,Z = polarization loss factor coordinates " 
o, = incoming wave unit polarization vector * 
tl, t2 = radiated power increments for @ and 6 
il = secondary field increments 


N,M = summation increments 
The first two inputs are physical dimensions’ obtained 
primarily through photographs, inputs with the “ are 
parameters that are either known or estimated, and the last 
inputs are only used to affect the processing time of the 
application. 

In order to calculate the far-field patterns of a 
parabolic reflector, the electric field at the aperture of the 
device must be known. Since the electric field patterns of 
interest are in the far-field, it is assumed that the circular 
aperture created by the reflector is located at the origin and 
is parallel to the z = 0 plane. The Cartesian components of 


the aperture electric field are expressed as [Ref. 14: pp. 


121-123; Ref. 15: pp. 29-31}: 
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Eax(0’,¢’) = [- Eg, (8’(r’) , 9’) cos(¢’) Ce 
eRe Ey, (8/ (r’) ,¢’)sin(¢’) esses etoey (V/m) 


Eay(r’, 9’) = [Ege (8 (r’) ,¢’) sin(¢’) ses 
4H Eae( 0 (me jy oe! yoos( oy jee ee estes (V/m) 


In (8.1) and (8.2), (E,,(r’,6’)) and (E,,(r’,0’)) are the 
aperture electric fields and k is the wavenumber (27/A). 
The relationships between (8’), (p), and (r’) needed to 


evaluate (8.1) and (8.2) are expressed as [{Ref. 4: p. 426]: 


Q/ = 2ean| | ee, (8.3) 


_ 4(foc)*+(r')? 


8.4 
4foc (m) ( 


p (r’) 

The electric field from the entire parabolic assembly is 
commonly referred to as the secondary pattern. Referring to 
Figure 8.1, the electric field in the -z half-plane is 
considered the secondary pattern. The secondary field of a 
parabolic reflector is calculated using the electric field 
integral equation solution and computing the component vector 


potentials (P,(6,¢)), (P,(8,¢)) as follows [Ref. 4: p. 426]: 


P,(8, o) -{"f fret Bar, ’) [e7kr’sin(8)(cos (61 cos () + sin($’) sin($))] 
aad (8.5) 


*« r/dr' dd’ (V/m) 
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2a a £ : weer , : , ; 
Ps (8 ; p) = J rs E.y (re ; d’) e jkr sin(@)(cos(¢ )cos(¢) + sin(¢ )sin(¢)) 
cae. (8.6) 


*«r/ dr’ dd’ (V/m) 


Double integrals in Mathcad take significantly longer to 
converge to an answer than do single integrals [Ref. 15: p. 
200}. To reduce the processing time of the applications, 
(8.5) and (8.6) were each separated into two single integrals. 
Therefore, for the applications, (8.5) and (8.6) are each 


computed as: 


P(8,o) = ['"A(8,0,6) db’ (V/m) (8.7) 
0 


where 


A(6,o, 6’) = f°" E(x’, 6’) @ dkr’sin (8) (cos (") cos () + sin(o”) sin(o)] -/ gy / (8.8) 
0 


Any contribution from feed system back lobes is ignored, 
and the component vector potentials calculated in (8.7) and 
(8.8) are applied to the subsequent formulas to compute the 


secondary field [Ref. 4: p. 383]: 


; ~Ikrg¢2 
By (0,) = 2X2 “1p (8, ) cos (b) + P, (8,6) sin (6) ] (v/m) (8.9) 
20 I ep 
. ake shrete _ . 
Ey (8,6) = =———_—_[-P, (0, 6) sin() + P, (8,6) cos ()]Jcos(8) (8.10) 


2U Ll eeo 


89 


The radiation intensity (U(8,¢)), radiated power (P,,4), 
and directivity (D,) of the secondary field pattern are 
expressed using the intrinsic impedance of free space (n, = 


1207) as follows [{Ref. 3: pp. 27-30]: 


2 
U(8,) = S12" 155 (8, 0) |? + |E, (8,6) |?]  (W/solid angle) (8-11) 
npn/ ; 
Pie f ih *u(8,o)sin(8) dbdb (Ww) (8.12) 
- = 27 Unax (dimensionless) (8.13) 
Prad 
In (8.13), (U,,,) is the maximum value of radiation intensity, 


obtained by creating (U,,,) aS a matrix and finding the maximum 
value of the matrix. 

The computer processing time required to evaluate (8.12) 
is prohibitively long. Thus, the parabolic reflector Mathcad 
applications use the following numerical approximation to 
evaluate radiated power [Ref. 3: pp. 37-42]: 


P.og = ETeay aps u(8,,,) sin(8,) (W) (8.14) 
rad N 2M . A’ ¥B A 


A=1 


In (8.14), (N and M) are user selected integers. This 
approximation for P,,, reduces the processing time of the 
applications substantially. As N and M increase, the results 


obtained with (8.14) converge to that obtained with (8.12). 
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When N = M = 3, the processing time for only P,,, of the 
helical feed reflector option is approximately four hours on 
a 33 MHz, 386 machine, and when N and M is increased to 5, the 
time to calculate P,,, for the same application is ten hours. 
For the applications, we chose N = M = 10 to evaluate P,,, with 
reasonable accuracy. With this choice, it takes approximately 
one day for the application to run, but this is much less than 
the thirty-six hours it takes if (8.12) is used. 

Parabolic reflectors are members of the aperture antenna 
family. The gain of an aperture antenna is the product of its 
aperture efficiency (e,,) and directivity. Therefore, the gain 
of a parabolic reflector is written as [Ref. 6: p. 575): 


G (dimensionless) (8.15) 


u 
o 
& 
% 
° 


Aperture efficiency is the product of several separate 
terms, including spillover efficiency (¢€,), taper efficiency 
(e€,.), phase (or random surface error) efficiency (e€,), 
polarization efficiency (¢€,), blockage efficiency (€,), spar 
efficiency (e€,,), and ohmic efficiency (€ omic) - 

Spillover efficiency (e,) is the fraction of power 
radiated by the feed system that is intercepted and reflected 
by the parabolic dish. Spillover efficiency is the most 
important efficiency term for any feed system and is defined 


as [Ref. 6: pp. 583]: 
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[?" [eu ¢(6") sin (6") ab! ae! 
ie) 0 


€ = 


_ (dimensionless) (8.16) 


[f° 0¢(8") sin (8!) a0’ ao! 
Q [8 


In (8.16), (U,;(8’)) is the radiation intensity function of the 
feed antenna and (8,) is the value of (8.3) when (r’) equals 
the radius of the mouth of the reflector (a,,,). 

Taper efficiency (e,) accounts for the lack of amplitude 
uniformity of the feed pattern on the surface of the 
reflector. Taper efficiency, which is a dimensionless 


quantity, is expressed as [Ref. 3: pp. 626-627]: 


en a 8) a0" ae 
je (8.17) 


2n a : / / ! 
[ ib ey dd 


rad 


€ _-2cor( 2) 
2 


Phase efficiency (e€,) is a function of the power loss 
that occurs if the field at the mouth of the dish is not in 
phase at every point in the aperture. Phase efficiency is 
approximately given by [Ref. 4: pp. 434-435]: 


-4n6'10°a,., 


T (8.18) 


(dimensionless) 


The power that is lost due to cross-polarized electric 
fields in the aperture of the reflector determines 


polarization efficiency (e,). Polarization efficiency is 
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difficult to precisely calculate, but can be reasonably 
approximated as 0.98 [Ref, 4: p. 435]. 

Blockage (€,) and spar (€,,) efficiencies are determined 
by the presence of feed structures, support struts, and 
possible signal processing equipment in front of the 
reflector’s aperture. Precise computation of blockage and 
spar efficiencies requires an extensive method of moments 
simulation and is not attempted by the Mathcad parabolic 
reflector applications. However, if the radius (a;) of the 
blocking structure at the focal point and the number of spars 
(N) that are (A/2) thick is known, blockage efficiency and 
spar efficiency can be estimated using Tables 10.1 and 10.2, 


respectively [Ref. 4: p.436]. 


TABLE 8.1 Blockage Efficiency 
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Ohmic efficiency (€ m-) for a reflector is usually very 
high but difficult to accurately calculate. The Mathcad 
parabolic reflector applications assume ohmic efficiency is 
0.98. 

With all reflector efficiency terms known or estimated, 
the aperture efficiency (e,,), effective aperture (A,), and 


isotropic radiated power (EIRP) are given by [Ref. 3: pp. 623- 


630]: 
Os5e Eel pC Cyt, © sc © sunie (dimensionless) (8.19) 
Ag= €ap™(Azer) (m?) (8.20) 
EIRP = P,agD,  (W) (8.21) 


The magnitude of the feed current (I,) for the source 
antenna cannot be estimated based on measured geometry. 
However, if the feed current is normalized to 1 amp, radiation 


resistance (R,) for the entire reflector structure is: 


R, = ——22—(Q) (8.22) 


In (8.22), (P,,4) is the radiated power from the aperture of 
the reflector, not of the feed antenna. 
The Mathcad parabolic reflector applications use the 


radiation resistance computed for the parabolic reflector 
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antenna system to determine effective height (h,,) by [Ref. 6: 


p.- 42]: 


h.=2 r‘*e (m) (8.23) 


The bandwidth of a parabolic reflector is assumed to be 
the same as its feed antenna. 

The polarization loss factor (PLF) of the parabolic 
reflector for a user selected point in the secondary far-field 
and a specific wave unit polarization vector (o,) is found 
using the electric field components of (8.9) and (8.10) as 


follows [Ref. 3: 51-53]: 


PLF = |G,°63|7, (dimensionless) (8.24) 


where the unit polarization vector of the antenna is: 


aE, + ay Ey + a, 


VIE(X, yy, Z) | 


OAxX,Yi2Z) = (dimensionless) (8.25) 


and the Cartesian components of the electric field are [Ref. 


3: p. 776): 
E,= E,cos (8) cos() - E,sin(o) (V/m) (8.26) 
E,= Ecos (8) sin() + E,cos () (V/m) (8.27) 
E,= -E,sin(6) (V/m) (8.28) 
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All of the secondary field parameters calculated by the 
Mathcad parabolic reflector applications assume that the 
observation point is in the far-field. Thus, all of the 
subsequent formulas must hold for the application calculations 


to be valid [Ref. 4: pp. 24-25]: 


r21.6A (m) (8.29) 
r2z104,,, (m) (8.30) 

2 
rs Cree) (dimensionless) (8.31) 


Table 8.3 is a comparison of reflector directivity 
computed by the Mathcad applications to predictions made for 
an identical antenna with unity aperture efficiency [Ref. 3: 
p.- 634). As can be seen, the comparisons vary by 
approximately 8 dB. This can be attributed to the fact that 
Prag WAS approximated using (8.14) instead of (8.12), and also 
to the fact that the numerical tolerance variable (TOL in 
Mathcad’s math menu) for the application was increased to 0.5 
verses 0.001 in order to reduce processing time. This 
increase in TOL reduces both the accuracy of (8.14) and the 


processing time of the application. 
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TABLE 8.3 Parabolic Reflector Data Comparison 


ANTENNA PREDICTED CALCULATED 


PARAMETER DATA DATA 


DIRECTIVITY 


Figure 8.2 is a comparison of measured data fed by an 
optimum horn with that calculated by the Mathcad application 
fed by an equiangular spiral for a parabolic reflector (a,,; = 
0.320 m, f = 0.147 m) operating at 35 GHz [Ref. 14: p. 126]. 
As can be seen, the comparison is good. The Mathcad pattern 
of Figure 8.2 is not symmetric which can be attributed to the 
number of increments used to increase the speed of processing. 

The reflectors fed by pyramidal and conical horn antennas 
are not included in this report. The applications have been 
developed and are included on the diskette containing the 
software. The complexity of the software results in 
prohibitively long processing time for the two applications 
(i.e., the pyramidal horn fed reflector application was still 
calculating P,,4. after seven days when run on a 50 MHz, 486 


personal computer using Mathcad 4.0). 
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0.995 


-0.1 0.08 -0,06 -0.04 -0.02 0 0,02 0.04 0.06 


- Mathcad Pottern, ... Measured Pattern 


FIGURE 8.2 Comparison of Parabolic Reflector E-Plane 
Electric Field Patterns 
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IX. REMARKS AND CONCLUSIONS 


This thesis and the accompanying software applications 
are intended to provide the Naval Maritime Intelligence Center 
(NAVMARINTCEN) with software for the analysis of various 
antenna types where information regarding the antenna is 
restricted primarily to the antenna type and the physical 
dimensions of the antenna. The goal of this research was to 
provide user friendly software that an engineer familiar with 
basic antenna types can easily use and interpret. Although 
the software was developed for that purpose, the complexity of 
formulas that apply to many of the antenna types of necessity 
reduced the simplicity of the applications. 

The software applications developed are written in 
Mathcad 3.1 and are compatible with any IBM personal computer 
that supports Mathcad 3.1 or Mathcad 4.0 for Windows. Many of 
the applications require extensive processing time. In order 
to reduce processing time, which results from the complexity 
of various integrals, approximate numerical techniques were 
implemented to evaluate the integrals and the numerical 
tolerance of the more complicated applications was reduced. 
Even with the processing time reduced, some programs require 
up to five days to complete when run with Mathcad 4.0 ona 50 
MHz, 486 personal computer. Applications that did not exceed 


a day to complete were run on a 33 MHz, 386 personal computer 
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with Mathcad 4.0. All applications were run using Mathcad 4.0 
Since Mathcad 4.0 iS approximately twice as fast as Mathcad 
3.1. As previously mentioned, the applications were written 
in Mathcad 3.1 due to sponsor requirements. Mathcad 3.1 is 
compatible with Mathcad 4.0, but Mathcad 4.0 is not compatible 
with Mathcad 3.1. 

For this thesis, a number of the antennas could not be 
adequately researched through the use of texts. In some 
cases, such as the caged dipole, very little information on 
the antenna could be found in either texts or in professional 
journals. For the caged dipole, the antenna is modeled as a 
circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles to take into account the number of 
conductors that are present in the caged dipole. This model 
provides good results when compared to predicted (the caged 
dipole’s performance is similar to that of a thick cylindrical 
dipole) and ELNEC results. 

It was found that Mathcad is not the ideal tool to 
analyze antennas due to its long processing time. An 
additional drawback of Mathcad is that when a function is in 
the beginning of an application and the same function is 
placed later in the application, Mathcad does not store the 
values of the function, but re-evaluates the function. This 
re-evaluation increases processing time significantly, 
especially when the function itself takes a very long time 


(hours vice minutes) to reach a numerical solution. 
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Nevertheless, the Mathcad applications developed in this 
research allow the user to analyze several antenna types 
(i.e., arrays with isotropic point sources, parabolic 
reflectors with the various feed antennas) that could not 
easily be evaluated using some antenna numerical analysis 
programs such as NEC, ELNEC, or MININEC. For the other 
antenna types researched, the antenna numerical analysis 
programs mentioned would be easier to setup, run, and provide 
parameters closer to those of a real antenna. The 
disadvantage of this approach is that it requires an ability 
to evaluate antennas using a program such as ELNEC, while the 
Mathcad applications do not require any advanced programming 


knowledge on the part of the user. 
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LINEAR ARRAY 
MATHCAD SOFTWARE ARRAY LN.MCD 


This linear array application is based on N isotropic point sources 
uniformly spaced and aligned on the z-axis. Mutual interference between 
adjacent sources are ignored when calculating antenna parameters. Since 
isotropic point sources are used as the radiating elements, polarization 
of the antenna can not be determined. The gain can only be idealized 
when it is assumed that the antenna input resistance is equivalent to 
the radiation resistance. Bandwith calculations are not carried out as 
a result of the narrowband characteristics of arrays, and in large part 
to the dependance inter-element spacing (d) has with the wavelength (A). 
The following limits should be adhered to when analyzing Broadside, 
End-Fire, and Phased (scanning) antennas: 


0 << 2 (Scanning Antennas) 
0) = 1/2 (Broadside Antennas) 
8) = O or n (Ordinary End-Fire Antennas) 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
Parameters do not contain subscripts (i.e., AF, will be written AFn).) 


The linear array antenna Mathcad application will compute the following 
parameters (Items with * indicate parameters that can be computed but 
their results are idealized as the result of using isotropic point 
sources as the radiating elements): 


} = Wavelength 

Imin = Minimum Distance to the Far-Field 
B = Phase Difference Between Adjacent Elements 
y = Array Factor Phase Shift 

AF, = Normalized Array Factor 

U(8) = Radiation Intensity 

Praq = Radiated Power 

Do = ODirectivity 

EIRP = Effective Isotropic Radiated Power 
Rr, = Radiation Resistance 

Rin = Antenna Input Resistance * 

Tl = Voltage Reflection Coefficient * 

tf = Antenna Efficiency * 

G = Gain * 

Rem = Maximum Effective Aperture * 

hem = Maximum Effective Height * 
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The following known or estimated data must be entered: 


N = Number of Isotropic Radiating Elements 
f = Frequency of Interest 
ad = Distance Between Adjacent Elements 
64 = Direction of Main Lobe 
In = Antenna Feed Current 
Zo = Characteristic Feed Impedence 
Enter input data here 
N:= 10 (elements) d:= 0.25 (meters) 
f:=3-10° (Hz) Go := 0 (radians) 
Io = 1.0 (A) Zo := 75.0 (Q) 


Define constants and calculate wavelength: 


¢ := 2.9979-10° (meters/sec) no := 120-1 (Q) 
at 2 or = 4 
ASH (meters/cycle) k :S ae (m7 +) 
f rs 
4 =0.999 (meters/cycle) ecd ‘= | (dimensionless) 
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Calculate linear array parameters: 


Define angular offset @: 


8:5 (o- 104 (2) F (2-2- 104) (radians) 


Minimum Distance to Far-Field rnin: 


r= 1.6-2 (meters) 
r :=5-N-d (meters) 
2 
r 22 (Nd) (meters) 
2 
a 
rmin := max(r) (meters) 
rmin = 12.509 (meters) 


Array Factor AF nt 


B :=-k-d-cos( 60) (radians) 


y(@) = k-d-cos(@) + B (radians) 


sin yo 
Agate) S28 


sn wo) 
2 


(dimensionless) 
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Radiation Intensity U: 


U(6) *= (AFn(8))° (W/solid angle) 


Radiated Power Prad: 


i 
Prad ‘= 2-1: U(8)-sin(6) dé (W) 
“0 


Prad = 1.256 (W) 


Directivity Do: 


Umax .= 1 (W/solid angle) 

Do := 4-n Umax (dimensionless) 
Prad 

Do = 10.007 (dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad-Do (W) 


EIRP = 12.566 (W) 
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Radiation Resistance R;: 


Rris zene 
([Iol ) 
Rr = 2.512 


Input Resistance Rip: 


Rin := Rr 


Rin = 2.512 


Voltage Relection Coeffecient [: 


rez Rin - Zo 
Rin + Zo 


lr =-0.935 


Reflection Efficiency ¢e: 


a zeed[1- (Irl)?] 


et = 0.125 
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(Q) 


(Q) 


(Q) 


(Q) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain G: 


G.=et-Do (dimensionless) 


G = 1.255 (dimensionless) 


Maximum Effective Aperture Rem: 


2 
Aem ‘= Ga 


4-7 


Aem = 0.1 


Maximum Effective Height Nem: 


hem := 2. Rr-Aem 
no 


hem = 0.052 
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GdB .= 10-log(G) 


GdB = 0.986 


(dB) 


(daB) 


(m) 


NORMALIZED ARRAY FACTOR RECTANGULAR PLOT 


Number of radiating elements: 
N =10 
Direction of Main Lobe (radians): 


60 = 0 


0.75 


Larni @)| 0.5 


0 0.628 1.257 1.885 2.513 3.142 3.77 4.398 5.027 5.655 6.283 
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NORMALIZED ARRAY FACTOR POLAR PLOT 


Number of Elements: =10 
Direction of Main Beam {radians): Go = 0 

n ‘ n 

X(0) i= [AFn(6)] -cosl 6 + — ¥(0) = [AFn(@)|-sino+ © 


Sane Sees 


Fe a a a 


sy 


' 

’ 
~~ 

hg = 


‘ 
‘ 
‘ 
' 
‘ 

rc 
‘ 
‘ 
‘ 
‘ 
¢ ‘ 
' ‘ 
U ) 
“re Soe 
‘ ‘ 
' 
‘ 
‘ 
‘ 
‘ 
‘ 
L 
’ 
‘ 
‘ 
‘ 
‘ 
‘ 
‘ 
‘ 
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PLANAR ARRAY 
MATHCAD SOFTWARE ARRAY _PL.MCD 


This planar array application is based on M and N isotropic point 
sources uniformly spaced and aligned on the x-axis and y-axis 
respectively. Mutual interference between adjacent elements are 
neglected when calculating antenna parameters. Since isotropic point 
sources are used as the radiating elements, polarization of the antenna 
can not be determined. The gain can only be idealized when it is 
assumed that the antenna input resistance is equivalent to the radiation 
resistance and the total antenna efficiency is calculated for a ideal 
lossless antenna. 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (i.e., AF, will be written AFn).) 


The planar array antenna Mathcad application will compute the following’ 


parameters (Items with * indicate parameters that can be computed but 
their results are idealized as the result of using isotropic point 


sources as the radiating elements): 


A = Wavelength 

Ly = Antenna Dimension in x-direction 

Ly = Antenna Dimension in y-direction 

Imin = Minimum Distance to the Far-Field 

By = Progressive Phase Shift Between Adjacent Elements in x-direction 
By = Progressive Phase Shift Between Adjacent Elements in y-direction 
Yyx = Array Factor Phase Shift in x-direction 

Yy = Array Factor Phase Shift in x-direction 

AF, (0,6) = Normalized Array Factor 

U(06) = Radiation Intensity 

Praq = Radiated Power 

Do = Directivity 

EIRP = Effective Isotropic Radiated Power 

Rr = Radiation Resistance 

Rin = Antenna Input Resistance * 

[T = Voltage Reflection Coefficient * 

€& = Antenna Efficiency * 

G = Gain * 

Rem = Maximum Effective Aperture * 

Nem = Maximum Effective Height * 
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The following known or estimated data must be entered: 


M = Number of Isotropic Radiating Elements in x-direction 
N = Number of Isotropic Radiating Elements in y-direction 
f = Frequency of Interest 
d, = Distance Between Adjacent Elements in x-direction 
dy = Distance Between Adjacent Elements in y-direction 
8, = Direction of Main Beam at @ = @o0 
$, = Direction of Main Beam at $ = $0 
Ig = Antenna Feed Current 
Zo = Characteristic Feed Impedence 
} = Number of Increments in Degrees for Far-Field Radiation Pattern 
Enter input data here 
M :=5 (elements) dx :=0.5 (meters) 
N =5 (elements) dy =0.5 (meters) 
. 8 1 ; 
f:=3-10 (Hz) Qo ‘= — (radians) 
6 
n 
lo = 1.0 (A) go =— (radians) 
Zo := 75.0 (Q) j:=180 (degrees) 
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Define constants and calculate wavelength: 


c = 2.9979-10° (meters/sec) no = 120-1 
_€ _2n 
Vise (meters/cycle) k= — 
f d 
4 = 0.999 (meters/cycle) ecd ‘= ] 


Calculate linear array parameters: 


Define angular offset @ andg: 


9 =(o- 104 { 1044 = one 104 
j 


b:=(0- 104 { 1044 2 Qa 104 


Minimum Distance to Far-Field Imin: 


r= 162 (meters) 


122 5-4(M-dx)?+ (N- dy)? (meters) 


3° (meters) 


A 


rmin ‘= max(r) (meters) rmin = 25.018 
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(Q) 


(m71) 


(dimensionless) 


(radians) 


(radians) 


(meters) 


Array Factor AFyt 


Bx :=- k-dx-sin(80)-cos( $0) 
By =-k-dy-sin(@0)- sin( $0) 
wx(6, o) ‘=k-dx-sin(@)-cos( >) + Bx 


yy( 8,9) "=k-dy- sin(@)-sin(>) + By 


sn wc. snl 00,0 
AFn(6,) ‘= _\2 : ee 
M-sin [ecw N-sin een) 

2 2 
Radiation Intensity U: 


U(6, 6) = AFn(6, $)- AFn(6, 6) 


Radiated Power Prad: 


2-7 = 
2 
Prad ‘= U(6,)-sin(6) dé do 
‘QO 0 
Prad = 0.206 


(radians) 


(radians) 


(radians) 


(radians) 


(dimensionless) 


(W/solid angle) 


(W) 


(W) 
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Directivity Do: 


Umax := 1 (W/solid angle) 

Do:= 4x Umax (dimensionless) 
Prad 

Do =61.09 (dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP = Prad-Do (W) 


EIRP = 12.566 (W) 


Radiation Resistance Ryi 


Rrcs aoe (Q) 
(IIo ) 

Rr = 0.411 (Q) 

Input Resistance Rint 

Rin := Rr (Q) 

Rin = 0.411 (Q) 
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Voltage Relection Coeffecient [: 


r 2 Rin- Zo (dimensionless) 
Rin + Zo 
lr =-0.989 (dimensionless) 
Reflection Efficiency e&: 
: 2 F : 
et ‘= ocd-| l- dr) (dimensionless) 
et = 0.022 (dimensionless) 
Gain G: 
G ‘=et:Do (dimensionless) GdB ‘= 10-log(G) (dB) 
G = 1.326 (dimensionless) GdB = 1.225 (aB) 
Maximum Effective Aperture Agm: 
Gv 
Aem != (m2) Aem =0.105 (m2) 
4-n 


Maximum Effective Height Rem? 


hem ‘= 2. Reacm (m) hem = 0.02! (m) 
| n0 
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TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 


Number of Elements along x-axis: M=5 
Number of Elements along y-axis: N=5 
Direction of Main Beam: 80 = 0.524 (radians) 
go = 0.785 (radians) 
X1(6) ‘= |AFn(6,0)] cos ine Y1(6) ‘= |AFn(0,0)|-sin]o- 2 
@ = 0° 


y 
“eo 
ra 
c 
= ' 
7 
: ’ 
’ 


. poets pee des ae aus 


= - 
~ 
.* 


~ 
sn 
‘N 


4 


0 
xref, X1(8) 
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TWO-DIMENSIONAL ANTENNA PATTERN: (y-2) plane 


Number of Elements along x-axis: M=5 
Number of Elements along y-axis: N=5 
Direction of Main Beam: Bo = 0.524 (radians) 
go =0.785 (radians) 
X2(6) = |AFn|6,—| | -cos0— — Y2(6) =|AFn|6,—}|-sinfo- — 
2 2 2 2 
§ = 0° 


ore beavete a eres. Se 


<8 
oo 

er 
. 


xref X2(9) 
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TWO-DIMENSIONAL ANTENNA PATTERN: $ = n/4 


Number of Elements along x-axis: M=5 
Number of Elements along y-axis: N=5 
Direction of Main Beam: Qo = 0.524 (radians) 
go =0.785 (radians) 
X3(0) := |AFn[0,—] |-cos]o- — Y3(0) := | AFn|0,—| | -sin[o- = 
4 2 4 2 
@ = 0° 


eS Se SSS See ? we ee eb - - -- eee em me mee eer err 


a ee 


a+? 


Se ee ae ee 


Se. 


= 0 1 
xref, X3(8) 
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TWO-DIMENSIONAL ANTENNA PATTERN: ¢ = by, 


Number of elements along x-axis: M=5 
Number of Elements along y-axis: N=5 
Direction of Main Beam: Go = 0.524 (radians) 
do =0.785 (radians) 
X4(6) := |AFn(0, $0)| -cos|@— ; Y4(0) := | AFn(6,40)| -sin| 6 - = 
6 = 0° 


oe ee ee 


Pad 


oit'=. 


- 


’ 
‘. 
~~ 


-] 0 1 
xref, X4(0) 
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CIRCULAR ARRAY 
MATHCAD SOFTWARE ARRAY _CI.MCD 


This circular array application is based on N isotropic point sources 
uniformly spaced on the x-y plane along a circular ring of constant 
radius. Mutual interference between adjacent elements are neglected 
when calculating antenna parameters. Since isotropic point sources are 
used as the radiating elements, polarization of the antenna can not be 
determined. The gain can only be idealized when it is assumed that the 
antenna input resistance is equivalent to the radiation resistance and 
the total antenna efficiency is calculated for a ideal lossless antenna. 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (i.e., AF, will be written AFn).) 


The circular array antenna Mathcad application will compute thefollowing 
parameters (Items with * indicate parameters that can becomputed but 
their results are idealized as the result of usingisotropic point 
sources as the radiating elements): 


A = Wavelength 

Imin = Minimum Distance to the Far-Field 
a, = Phase Excitation of the nth Element 
AF (8,6) = Normalized Array Factor 

U(0,o) = Radiation Intensity 

Praq = Radiated Power 

Do = Directivity 

EIRP = Effective Isotropic Radiated Power 
Ry = Radiation Resistance 

Rin = Antenna Input Resistance * 

[ = Voltage Reflection Coefficient * 

f& = Antenna Efficiency * 

G = Gain * 

Rem = Maximum Effective Aperture * 

Nem = Maximum Effective Height * 
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The following known or estimated data must be entered: 


N Number of Isotropic Radiating Elements 
f = Frequency of Interest 
a = Radius of circle 
8) = Direction of Main Beam at @ = @o 
¢g = Direction of Main Beam at 6 = go 
Ip, = Antenna Feed Current 
Zo += Characteristic Feed Impedence 
} = Number of Increments in Degrees for Far-Field Radiation Pattern 
Enter input data here 

. _ 10 

N := 10 (elements) ais (meters) 
2-0 

f:=3-108 (Hz) Go := 0 (radians) 
Io = 1.0 (A) go := 0 (radians) 
Zo := 75.0 (Q) j:=360 (degrees) 
Define constants and calculate wavelength: 
¢ = 2.9979-108 (meters/sec) no :=120-n (Q) 

_c¢ een -1 
oe (meters/cycle) ki=—— (m7 *) 

r 
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Define constants and calculate wavelength: 


4 =0.999 (meters/cycle) ecd ‘= | 


Calculate linear array parameters: 


Define angular offset § and: 


e:=(0- 104), 10°44 2. (22-10) 
j 


6 =(0- io), ioe 1. (2-2- 10°4) 
j 


Minimum Distance to Far-Field Ymins 


Ty =1.6-2 (meters) T =5-(2-a) 
2 
rea) (meters) 
r 
rmin := max(r) (meters) rmin = 20.278 


Array Factor AF: 


n=l1.N (elements) ®@ =2-1— 


acs: k-a-sin(80)-cos{ 4o - o) 


122 


(dimensionless) 


(radians) 


(radians) 


(meters) 


(meters) 


(radians) 


(radians) 


Array Factor AF: 
AF(6,$) =< x - (k-a-sin( 0)-cos (9 - On) + a) 


n 


Radiation Intensity U: 


U(6,$) = AF(6,$)- AF(6, >) 


Radiated Power Praq: 


2:0 = 
Prad := U(6,$)-sin(@) d6 dd 
‘Oo +0 
Prad =0.536 


Directivity Do: 


Umax := 1 


_ 4-1 Umax 
Prad 


Do. 


Do = 23.427 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad-Do 


EIRP = 12.566 
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(dimensionless) 


(W/solid angle) 


(W/solid angle) 


(dimensionless) 


(dimensionless) 


Radiation Resistance R;: 


Rris aise (9) 
(lol ) 

Rr = 1.073 (Q) 
Input Resistance Rint 
Rin := Rr (Q) 
Rin = 1.073 (Q) 
Voltage Relection Coeffecient [: 
liz pu (dimensionless) 

Rin + Zo 
l =-0.972 (dimensionless) 
Reflection Efficiency et: 

2 ; . 

et = sed-| 1 (irl) | (dimensionless) 
et =0.056 (dimensionless) 
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Gain G: 


G.=et-Do (dimensionless) 


G = 1.303 (dimensionless) 


Maximum Effective Aperture Remi 


G12 


4-n 


Aem ‘= (m2 ) 


Maximum Effective Height Dem: 


hem .= 2: Rr-Aem (m) 
70 
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GdB = 10-log(G) 


GdB = 1.149 
Aem =0.104 
hem = 0.034 


(dB) 


(dB) 


(m2 ) 


(m) 


TWO-DIMENSIONAL ANTENNA PATTERN: (x-zZ) plane 


Number of Elements in x-y plane: N =10 
Direction of Main Beam: Go =0 (radians) 
go =0 (radians) 
: n ; 1 
X1(0) := |AF(6,0)| -cos|@— Y1(@) := |AF(0,0)|-sin ee 
6 = 0° 


¢ 
¢ 
é 
‘ 
Ly 
r 
i] 
‘ 
i] 
‘ 


-} 0 1 
xref .X1(8) 
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TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 
N =10 


Number of Elements in x-y plane: 


Direction of Main Beam: 


X2(6) = |AF es 


-cos]6- su 
2 


60 =0 


do =0 


Y2(6) =|AF 05 


(radians) 


(radians) 


-sin{@- — 
2 


MS ere ao ge teta rare apse Sei 
4 


¢ 


U 
4 
a 
4 


xref. X2(6) 


127 


TWO-DIMENSIONAL ANTENNA PATTERN: $ = 2 


= a/4 
Number of Elements in x-y plane: N =10 
Direction of Main Beam: 60 =0 (radians) 
go =0 (radians) 
X3(0) := |AFl0,—|| cosla- — Y3(0) := |AF[0,~]{-sin]o- — 
4 2 4 2 
@ = 0° 


¢ 
c 
' 
' 


a 
‘ 
‘ 
‘ 
L] 
¢ 
~ 
? 
t 
? 
4 
? 
‘ 
a 


= 
3 


0 
xref» X3( 8) 
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TWO-DIMENSIONAL ANTENNA PATTERN: $ = v5 
Number of Elements in x-y plane: N =10 
Direction of Main Beam: 80 =0 (radians) 
go =0 (radians) 
fog ' 1 
X4(6) *= |AF(0,4¢0)| -cos| 6— - Y4(0) := | AF(0,40)| -sin}@— : 
\ 
6 = 0° 


‘ 
i] 
‘ 
~~ 
? 
L] 
td 
’ 
é 
~ 


4 


xref» X4( @) 
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FOLDED DIPOLE 
MATHCAD SOFTWARE FOLDED.MCD 


To provide impedance matching to coaxial input cables and maintain good 
directional pattern characteristics, a folded dipole is used instead of 
a Single dipole element. The geometry of a folded dipole is such that 
the spacing between dipole elements (d) is much less than the wavelength 
(4) and the length of an element (L). 


This Mathcad application calculates the parameters of a folded dipole 
aligned vertically on the z-axis and a folded dipole positioned 
horizontally above the earth. For the calculations, the antenna is 
assumed to behave as a cylindrical dipole with equivalent radius (a,). 


{Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (i.e., Rj, will be written Rin).) 


The folded dipole antenna Mathcad application will compute the following 
parameters (* indicates antenna parameters that are calculated for a 


folded dipole in free space and positioned horizontally over the earth): 


i} = Wavelength 

ae = Equivalent Radius 

Rr = Radiation Resistance 

Rq = Input Resistance of Linear Dipole 

Xr, = Reactance of Linear Dipole 

Xq = Input Reactance of Linear Dipole 

Zq = Input Impedance of Linear Dipole 

Zo = Characteristic Impedance of Two-wire Transmission Line 
Ze = Characteristic Impedance Transform of Folded Dipole 
Zin = Characteristic Impedance of Folded Dipole 
Zo = Maximum Current of Folded Dipole 

Imin = Minimum Distance to the Far-Field 

E(8) = Electric Field * 

E,(8) = Normalized Electric Field * 

U(8) = Radiation Intensity * 

Umax = Maximum Radiation Intensity * 

Praq = Radiated Power * 

Do = ODirectivity * 

EIRP = Effective Isotropic Radiated Power * 

Rin = Antenna Input Resistance * 

[. = Voltage Reflection Coefficient * 
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f& = Antenna Efficiency * 


= Antenna Polarization Vector * 


~ 
ia 
ry 

Il 


Polarization Loss Factor * 
G = Gain * 


Rem = Maximum Effective Aperture * 

hem = Maximum Effective Height * 

Erp = Relative Complex Effective Dielectric Constant of Ground 
Ty = Vertical Reflection Coefficient of Ground Plane 

Ty, = Horizontal Reflection Coefficient of Ground Plane 


The following known or estimated data must be entered: 


f = Frequency of Interest 

a = Radius of Conductor 

b = Radius of Conductor 

ad = Distance Between Conductors 

L = Length of Antenna 

h = Height of Antenna Above Earth 


Eoqg = Conduction/Dielectric Efficiency 

fr = Relative Dielectic Constant of Ground 

6 = Conductivity of Ground 

r¢¢ = Far Field Distance 

1 = Number of Increments in Degrees for Far Field Radiation Pattern 
g 

Ow = Incoming Wave Electric Field Unit Vector For Free Space Antenna 

Owh = %INcoming Wave Electric Field Unit Vector For Horizontal Antenna 

Ooh = Antenna Polarization Direction For Horizontal Antenna 
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ow [= }|0 

0 
f:=3-10° 
a := 0.0005 
bi=a 
d:= 12.5-a 
L:=0.5 
h:=.5 


Enter input data here 


(dimensionless) 


(Hz) 


(meters) 


(meters) 


(meters) 


(meters) 


(meters) 


owh °= 


rff:= 1-10° 


ecd ‘= ] 


Q 
m 
a 


(dimensionless) 


(radians) 


(meters) 


(radians) 


(dimensionless) 


(siemens/meter) 


(dimensionless) 


Define constants: 


oa -=|0 (dimensionless) 
0 
¢ := 2.9979-10° (meters/sec) no = 120-1 
c 2-1 
A=- (meters/cycle) ki 
ry 
4 =0.999 (meters/cycle) y ‘= 0.57721 
«lias _ 
£0 -=——- 10 (E/m) zZ-=0..21 
36-1 


Define cosine [Ci(x)] and sine [Si(x)] integrals: 


n:=1..50 


(- 1 ae 1 yan 


Ci(x) -= 0.57721 + In(x) - Si(x) = 
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(Q) 


(m71) 


(dimensionless) 


(radians) 


(radians) 


(- i= fan 1) 
(2-n)-(2-n)! (2-n- 1)-(2-n= 1)! 
n n 


FOLDED DIPOLE ALIGNED ON THE Z-AXIS IN FREE 


Eguivalent radius ag: 


a’-In(a) + b*In(b) + 2-a:b-In(d) 


ae = exp 
(a+ b) 


ae = 0.002 


Electric field E(Q): 


k-L k-L 
ae cos ao) = cos 
-j ck-rff 2 2 


E(8) =j -jo<——— 
2-arff sin(9) 
En(@) -- E(8) 


Minimum Distance to Far-Field Imin: 


ry i= 16:2 (meters) r 
2 
r, -2C) (meters) 
ry 
rmin ‘= max(r) (meters) 
rmin =2.5 (meters) 
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SPACE: 
(meters) 
(meters) 
(V/m) 
(V/m) 
=5-L (meters) 


Radiation Intensity U(0): 
ue) := (E891)? (W/solid angle) 
2-70 


Radiated Power Prsaq: 


1 
Prad := 2-1- U(8)-sin( 8) dé (W) 
‘0 


Prad = 36.64 (W) 


Directivity Do: 


uinax:= uf (W/solid angle) 

2 

Do '= amex (dimensionless) 
Prad 

Do = 1.641 (dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad-Do (W) 


EIRP = 60.132 (W) 
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Radiation Resistance Ry: 


Rr ial 4+ In(keL) — Ci(keL) + t-sin(k-L)-(Si(2- KL) - 2-Si(KL)) .. 


+ Leo rt nf + Ci(2:k-L)- 2-Ci(k-L) 


Rr = 73.281 


Input Radiation Resistance of linear dipole Ry: 


Rd ‘= Rr 
fo() 
sin| —— 
2 
Rd =73.281 


Reactance of linear dipole Xri 


Xr :2.[2-Si(k-L) 4 cos(k-L):(2:Si(K:-L) - Si(2-KL)) .. 
4-n 2-k-ae | 


+} - 1-sin(k-L)-]2-Ci(k-L) - Ci(2-k-L) - Ci 


Xr = 43.142 


Input reactance of linear dipole Xqi 


ie eats 
( | (A) 
sin|--— 
2 
Xd = 43.142 
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(Q) 


(Q) 


(Q) 


(Q) 


(Q) 


(Q) 


(Q) 


(Q) 


Input impedance of linear dipole 2q: 


Zd =Rd+j -Xd (Q) 


Zd = 73.281 +43.1421 (Q) 


Characteristic impedance of two-wire transmission line 2): 


Zo := 2°. acosh 7 (2) 
R 2-afa-b 


Zo = 302.312 (02) 


Characteristic impedance transform for folded dipole Ze 


Zt [= ) Zoran(K4| (Q) 
Zt =-2.747 10°; (Q) 


Characteristic impedance of folded dipole 2jy,: 


in 442d (Q) 
2:-Zd + Zt 
Zin = 293.307 +172.4651 (Q) 
Input Resistance Rin: 
Rin = Re(Zin) (Q) Rin = 293.307 (Q) 
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Voltage Relection Coeffecient [: 


ae Rin - Zo 
Rin + Zo 


(dimensionless) 


Reflection Efficiency e&: 


at secd(1- CIrl)’ (dimensionless) 


Gain G: 
G:=et-Do (dimensionless) 
G =1.641 (dimensionless) 


Polarization loss factor PLF: 


—| 2 
PLF = (low-oal) (dimensionless) 


Maximum Effective Aperture Rem? 


Gar 


4-n 


Aem ‘= -PLF (m2) 


Maximum Effective Height Dem 


hance eee. iat 
yo 


lr =-0.015 


et =] 


GdB := 10-log(G) 


GdB = 2.15 
PLF =1 

Aem = 0.13 
hem = 0.318 
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(dimensionless) 


(dimensionless) 


(dB) 


(dB) 


(dimensionless) 


(m) 


Far Field radiation pattern of vertical folded dipole in free space. 


X(0) := |En(@)|-cos 042 Y(0) := [En(0)| -sinfo+ — 


wv oe? 


Frew 


t 
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HORIZONTALLY POSITIONED FOLDED DIPOLE OVER THE EARTH: 


Define angular offsets 01 and 02: 


61 el pea (radians) @2 eee, (radians) 


Relative complex effective dielectric constant Erpi 


erp ‘= er—j - ; (dimensionless) 
2-1-f-€0 
¢ component of E-field Eo: 
j ckerff cos Ein a) - co 5] 
CIR a een | ele eee Se 2 (V/m) 


2: a-rff : 2. 2 
1 = sin(61)-sin(0) 


kL. kL 
cos| —-sin[62_}-sin(O)| - cos|— 
ee | ana | | sl 


El) 40 ——— | (V/m) 


Qari fie sin(or}?.sincoy 


@ component of E-field EQ: 


ei ket [oo sin(61 )- sn 3] - cos} 


2arit 


E0( 61) :=j ee (V/m) 


2 
1 = sin(@1 Pin 
2 


kL k-L 
cos] ——-sin[{@2_}-sin - cos|— 
ou (s (02) ) a 


2-nrff 


E01, -=j -jo——— (V/m) 
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Vertical plane wave reflection coefficient [,,: 


roby -« SB-20001)= ep sing)” 


erp-cos( 01) + «ferp - sin(@1 i 


: erp-cos(62,) - = sin(e2,)” 
op sn(a) 


(dimensionless) 


Tvl ——— (dimensionless) 


erp-cos( 02} + | - sin 


Horizontal plane wave reflection coefficient [,: 


Th(01) “= cos(@1) — dese sin(61 


(dimensionless) 
seeieclese sin(01)° 
cos(62) - lap sin(e2,}* 
rh) —— (dimensionless) 
cos(02} + lp - sin(@2)° 
E-field (% component) with ground plane reflection Edgp? 
Eégp(@1) -=E¢(@1)-(ei ERE" 4 ryegty.e5 Khost 81)} (V/m) 
j -kK-h- 62. <j) -k-h- 
Eggp, := me oa ,, Thl-e- oN) | (v/m) 
Egn(61) = [Etep ODL (Vv/m) 


max(Edgp 1) 
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E-field (@ component) with ground plane reflection E8gp 


Edgp( 61 ) := E6( 61 )- (J ehreos( 1) _ ry91 -e7! Aehreos(81)) 


(J -kch-cos(02) _ Tvl Pe nee? 
z 


Edgpl, := en, 


|Eegp(61 | 
max( E@gp1) 


E@n(61 ) := 


Radiation Intensity Uh(@1l): 


uncer) =. [¢fe¢gpco1 yl)?+ CLEoep¢01 | 92] 
2-no 


Radiated Power Phrag: 


a 
p 

Phrad ‘= 2-n- Uh(61 )-sin(@1 ) d6l 
0 


Phrad = 70.415 


Directivity Dho: 


Ugp_max | ject (|Beepr, Edgp 1 


2-no 


+ (fer) | 


Uhmax = max(Ugp_max ) 
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(V/m) 


(V/m) 


(V/m) 


(W/solid angle) 


(W) 


(W) 


(W/solid angle) 


(W/solid angle) 


Directivity Dhgo: 


Dho = 4-n Uhmax 
Phrad 
Dho = 3.938 


Effective Isotropic Radiated Power EIRPh: 


EIRPh := Phrad: Dho 


EIRPh = 277.3 


Input Resistance RAi ni 


Rhin 122 rad 
sy 
sm|— 
2 
Rhin = 140.83 


Voltage Relection Coeffecient [H: 


rH:= Rhin- Zo 
Rhin+ Zo 
TH =-0.364 


Reflection Efficiency chy: 


cht secd1 -( IrH] | (dimensionless) 
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(dimensionless) 


(dimensionless) 


(W) 


(W) 


(Q) 


(Q) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain Gh: 


Gh := eht-Dho (dimensionless) GhdB := 10-log(Gh) (dB) 


Gh =3.415 (dimensionless) GhdB = 5.334 (dB) 


Polarization loss factor PLFv: 


-h 
Ehmax ‘= | E6n( ph) | (V/m) Lobes ‘= toor|=* + } 
0 
0 
gah := (dimensionless) 
Ehmax 
Lobes 
——_]42 
pLFh:: ( lewh-oahl) (dimensionless) 
PLFh =7.736 10° (dimensionless) 


Maximum Effective Aperture AVem: 


Gh-i? 


4-n 


Ahem := -PLFh (m2) Ahem = 2.09% 10 *(m2) 


Maximum Effective Height Aven: 


hhem := 2. (him: Ahem (m) hhem =0.018 (m) 
no 
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H-Plane Radiation Pattern (y-z plane) of Horizontal Folded Dipole Over 
The Earth: 


Y1(61) = |E6n(61)|-sin{61 +s X1(61) = |Eon(61)| -cos| 61 +s 


yref 


Y1(91) 


xref, X1( 01) 
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E-Plane Radiation Pattern (x-z piane) of Horizontal Folded Dipole Over 
The Earth: 


Y2(01) = |Een(¢1)|-sin{¢l - X2(01) ‘= |E6n(¢1 )] -cos{ 61 iS 


xref. X2( 01) 
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CAGED DIPOLE ANTENNA 
MATHCAD SOFTWARE CAGED DI.MCD 


This caged dipole application is modeled after a circular array with the 
radiating elements as dipoles. The dipoles are modeled after conductors 
with infinitely small radii and antenna length greater than 1/4. Mutual 
interference between adjacent elements are neglected when calculating 
antenna parameters. 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (i.e., AF, will be written AFn).) 


The caged dipole antenna Mathcad application will compute the following 
parameters: 


VERTICAL CAGED DIPOLE IN FREE SPACE: 


4 = Wavelength 

Imin = Minimum Distance to the Far-Field 
AF (86) = Normalized Array Factor 

E(06) = E-Field of Finite Length Dipole 
E+ (06) = Total E-Field 

Etr (86) = Total Normalized E-Field 
U(0o) = Radiation Intensity 

Pragqd = Radiated Power 

Dy. = -Directivity 

EIRP = Effective Isotropic Radiated Power 
Rr = Radiation Resistance 

Rin = Input Resistance 

T = Voltage Reflection Coefficient 

& = Antenna Effeciency 

G = Gain 

PLF = Polarization Loss Factor 

Rem = Maximum Effective Aperture 

Nem = Maximum Effective Height 

flow = Low Frequency 

fhigh = High Frequency 

BW = Bandwidth 


VERTICAL CAGED DIPOLE OVER EARTH: 


Erp = Relative Complex Effective Dielectric Constant 


AFv (81,6) = Normalized Array Factor 
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Ev (61,6) = E-Field of Finite Length Dipole 


T,,(@1) = Vertical Plane Wave Reflection Coefficient 
Ev; (61,6) = Total E-Field 

Ever (61,96) = Total Normalized E-Field 

Uv (61,¢) = Radiation Intensity 

PVraq = Radiated Power 

Dvg = Directivity 

EIRPv = Effective Isotropic Radiated Power 
Rv, = Radiation Resistance 

Rvin = Input Resistance 

Iv = Voltage Reflection Coefficient 

ev; = Antenna Efficiency 

Gv = Gain 

PLFvV = Polarization Loss Factor 

AVem = Maximum Effective Aperture 

hVem = Maximum Effective Height 


HORIZONTAL CAGED DIPOLE OVER EARTH: 


AFh(61,$1) = Normalized Array Factor 

Ty (81) = Horizontal Plane Wave Reflection Coefficient 
E6h (61) = 6 Component of E-Field 

E6h; (61,61) = 6 Component of Total E~Field 

E€h; p, (61,91) = § Component of Total Normalized E-Field 
Egh (61) = $ Component of E-Field 

Eghy (81,61) = $ Component of Total E-Field 

Eghey, (01,61) = § Component of Total Normalized E-Field 
Uh (@1,$1) = Radiation Intensity 

Phraq = Radiated Power 

Dhoh = Directivity 

EIRPh = Effective Isotropic Radiated Power 

Rhy = Radiation Resistance 

Rh;y, = Input Resistance 

TH = Voltage Reflection Coefficient 

cht = Antenna Efficiency 

Gh = Gain 

PLFhH = Polarization Loss Factor 

Ahem = Maximum Effective Aperture 

hhem = Maximum Effective Height 
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The following known or estimated data must be entered: 


N = Number of Dipoles 

a = Radius of Antenna 

L = Length of Antenna 

h = Antenna Height Above Ground 
f = Frequency of Interest 

I = Antenna Feed Current 

a 


° 
o = Input Impedance 
fr = Relative Dielectric Constant of Ground 
6 = Conductivity of Ground 
Ow = %Incoming Wave E-Field Unit Vector for Antenna (Free Space) 
SO, = Antenna Unit Polarization Vector (Free Space) 
OVy, = Incoming Wave E-Field Unit Vector for Vertical Antenna 
OV = Direction of Incoming Unit Vector (Vertical Antenna) 
oh, = Incoming Wave E-Field Unit Vector for Horizontal Antenna 
Shy = Direction of Incoming Unit Vector (Horizontal Antenna) 
} = Number of Increments in Degrees for Far Field Radiation Pattern 
Enter input data here: 
1 ] 
ow (=| 0 (dimensionless) oa := |0 (dimensionless) 
0 0 
] 
80-1 ; 
ovw =|0 (dimensionless) §vp =—— (dimensionless) 
80 
0 
0 
ohw :=|0 (dimensionless) Ohp :=0 (dimensionless) 
l 
N =8 (elements) Ilo = 1.0 (A) 
ais (meters) Zo :=75.0 (Q) 


co{— 
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L:=0.5 (meters) j:=90 


h.= 0.25 (meters) 6 = 10 


f:=3-10° (Hz) er = 15 


Define constants and calculate wavelength: 


¢ = 2.9979. 10° (meters/sec) no ‘= 120-1 
are 2s 
Lise (meters/cycle) k i= — 
f A 
In zie (A) ff = 10° 
N 
me. 9 
co '=$ ——- 10 (farads/m) ecd ‘= ] 
36-1 
A =0.999 (meters/cycle) z-20..2-j 


Calculate caged dipole parameters parameters: 


Define angular offset @: 


6°=0,—..2-1 (radians) $°=0,—.. 2-n 
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(increments) 


(siemens/m) 


(dimensionless) 


(02) 


(m71) 


(meters) 


(dimensionless) 


(increments) 


(radians) 


Minimum Distance to Far-Field ryjn: 


I) *= 1.6-2 (meters) 


rie 5-4L24 (2-0)? (meters) 


— 2+ 2-2) 


4 (meters) 
d 


rmin ‘= max(r) (meters) rmin = 2.795 


VERTICAL CAGED DIPOLE IN FREE SPACE: 


Array Factor AF(@): 


n:=1..N (elements) ®@ i= 2-n— 


j -k-asin( 8)-coslo— 
AF(8,9) =) ine a a 


n 


E-Field of finite dipole E(0): 


In-e 


eet cos{ Ecos - co 
E(8) =] -]}0————- ae ey 2 
2-nrff sin(6) 
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(meters) 


(radians) 


(dimensionless) 


(V/m) 


Total E-Field E; (@)& Normalized E-Field E;, (8): 


Etn(6, 6) --__Et(8,6) (V/m) 
Aves 
2 2 


Radiation Intensit 


pacvet ten rete s ey 


U(B) = 
U(8, 4) = chewe.9) y (W/solid angle) 
2:0 
Radiated Power Prag: 
2:1 v4 (W) 
Prad ‘= U(6, >)-sin(8) d6 do 
10 0 
Prad =0.441 (W) 
Directivity Do: 
Hf 4-n Umax . 
Umax :=Uj/—,0 (W/solid angle) Do :=————— (dimensionless) 
2 Prad 


Do = 1.543 


(dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad-Do (W) 


EIRP = 0.681 (W) 
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Radiation Resistance Ry: 


_ 2:Prad 
([In| >” 


Rr: (Q) Rr = 56.487 (Q) 


Input Resistance Rini 


Rin ees (Q) Rin = 56.487 (Q) 


Voltage Relection Coeffecient [: 


ree Rin - Zo 
Rin + Zo 


(dimensionless) r =-0.141 (dimensionless) 


Reflection Efficiency ¢: 


at ‘ecd-[1- (|r])?] (dimensionless) et = 0.98 (dimensionless) 
Gain G: 

G:=et-Do (dimensionless) GdB ‘= 10-log(G) (aB) 

G =1.513 (dimensionless) GdB = 1.798 (aB) 


Polarization Loss Factor PLF: 


—y,2 
PLE «= ([ow-oal) (dimensionless) PLF = 1 (dimensionless) 
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Maximum Effective Aperture Aen: 


Ga? 


4-n 


Aem ‘= 


-PLF (m2) Aem =0.12 (m2) 


Maximum Effective Height hen: 


hem := 2. Rr-Aem (m) hem = 0.268 (m) 
| no 


Define Frequency Relationship F(x): 


F(x) =-5.696 10° °-x+4 0.31481 (dimensionless) 


Frequency Range f)5y_and fhigh: 


flow :=f- re] -f (Hz) flow = 2.056 10° (Hz) 
fhigh := f+ re} 4 (Hz) thigh =3.944 10° (Hz) 


Bandwidth BW: 


BW =fhigh- flow (Hz) BW = 1.888"10° (Hz) 
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E~FIELD PATTERN: This far-field radiation pattern is for the vertical 
caged dipole in free space: 


Number of Dipole Conductors: N=8 
Radius of Antenna: a=0.125 (meters) 
X(6) -= |Etn(0, 4) -cos{0- — Y(6) -= |Etn(0,¢)] -sin{o- : 
2 
@ = 0° 


4 
I 


Se pero 
C 


trccccce 


i] 

4 

' t] 
~~ er regret ree 


-~s err re 
‘ 


~~ 


’ 
’ 
' 
7 


LoS 


VERTICAL CAGED DIPOLE OVER EARTH: 


Define angular offset $1 and §2: 


01 ee ea (radians) $2 siz (radians) 
9° 9°35 ae 

62 erat, (radians) 

BD Dy 


Relative complex effective dielectric constant Erpi 


ep =er-j - (dimensionless) 
2-1-f-0 
Array Factor AFv(@): 
j -k-a-sin( 61)-cos(o— ® 
AFv(61,$) =S ‘ine’ iy a(t ) (dimensionless) 
n 
j -k-a-sin(62.) - 2,—- 0 
AF vgp | =) Ine’ snl) eo) (dimensionless) 
n 
E-Field of finite dipole Ev(@1l): 
th : 
= eos 'gew 61) - cos{ kL 
Ine? kt MD \2 
Ew@l) =) -7o— |——————_————“— (V/m) 
2-nrff sin(61 ) 
hai cos oer 62.}'| - cos Me 
ee 2 2 
Evgp *=) 0 * | (V/m) 
2-nrff sin(62) 
Vertical plane wave reflection coefficient [,,(01): 
| P 2 
PCO jee Poe = wep sinh) (dimensionless) 


erp cos( 1) + ile - sin(6] y’ 
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_ ap-cos{ 62 ] - jap- sin(o2)° 


I'vh: ——— (dimensionless) 
2 

erp-cos{ 62 } + jerp- sin(02 ) 
Total E-Field Ev; (81) & Normalized E-Field Evty (01): 
Evt(61,) = Ev(1)-AFV(61,4)-(e! Kh-cos( 81) 4 Py g] ).e74 Kehvcos( 01)} (V/m) 

j -k-h-cos -j k-h-cos(62 
Evtgp, *= | Evep, Arve, {e col) -e ena ) | (v/m) 
Evtn( 61 ,o) _ Emon. gt (V/m) 
max(Evtgp) 


Radiation Intensity Uv(Ol): 


Uv(01,$) 2 ae | Evtcel ,6)])? 
2-no 


Uvl eae Eupl) 
7 2-no 


Radiated Power PVrad: 


2n]- 
v 
Pvrad ‘= Uv(61,$)-sin( 61) dél do 
0 0 
Pvrad = 0.483 


Di rectivity Dvo: 


Uvmax ‘= Uv( @vp,0) 


(W/solid angle) 


(W/solid angle) 


(W/solid angle) 
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Directivity Dv,5: 


4: 1 Uvmax 
Dvo ‘= 
Pvrad 
Dvo = 4.415 


Effective Isotropic Radiated Power EIRPv: 


EJRPv ‘= Pvrad:Dvo 


(W) 


Radiation Resistance Rvypi 


__ 2-Pvrad 
(In|)? 


Rvr 


(Q) 


Input Resistance RVini 


Rvr 
Rvin ‘= 


(Q) 


Voltage Relection Coeffecient [Vi 


rV cs Rvin- Zo 
Rvin + Zo 


(dimensionless) 


Reflection Efficiency tvy: 


evt :=0ed| 1 - ({rvl)? 


(dimensionless) 
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EIRPv = 2.133 


Rvr =61.857 


Rvin = 61.857 


rV =-0.096 


evt = 0.99] 


(dimensionless) 


(dimensionless) 


(W) 


(Q) 


(Q) 


(dimensionless) 


(dimensionless) 


Gain Gv: 


Gv =evt-Dvo (dimensionless) 


Gv = 4.374 (dimensionless) 


Polarization Loss Factor PLFv: 


Evmax ‘= Evtn(@vp,0) (V/m) 


Evmax 


lobesv 
ova:= 
0 


0 


(dimensionless) 


ees 11.2 
PLFv = A Sea ) (dimensionless) 


Maximum Effective Aperture Agm: 


2 
Avem ‘= ove -PLFv (m ) 
4-n 


Maximum Effective Height hVen: 


hoemtee. Rvr-Avem an) 
no 


GvdB =10-log(Gv) (dB) 


GvdB = 6.409 (dB) 


lobesv ‘= floor + + } 
ry 


PLFv = 0.986 (dimensionless) 
Avem = 0.343 (m2 ) 
hvem = 0.474 (m) 


159 


E-FIELD PATTERN: This far-field radiation pattern is for the vertical 
caged dipole over earth: 


Number of Dipole Conductors: N=8 (elements) 
Radius of Antenna: a=0.125 (meters) 
Height of Antenna Above Earth: h=0.25 (meters) 
X1(61) := |Evin(61 ,4)|-cos{@l += Y1(01) := [Evin(61,6)]-sinfo1 + — 
2, 2 
8 = 0° 


‘ 

' 

' 

‘ 

‘ 

‘ 

‘ 

' 
cr 
‘ 

' 

' 

' 

‘ 

' 

' 
od 


xref» X1(81) 
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HORIZONTAL CAGED DIPOLE OVER EARTH: 


Define angular offset $1 and 92: 


(1 EO2 24 (radians) 
j 


Array Factor AFh(@1,$1): 


n‘=1.N (elements) © '2 2-1 


jcka (sinc $1 )-sin( 61 )-sin (0) +cos( 61 )-cos(@,}) 
€ 


AFh(O1,§1) = > In: 


n 


AFhI, != Stine’ -k-a-(sin(42) -sin(02) -sin{0,) + cos (62) -cos (0,)) 


@ component of E-Field of finite dipole E@h(@1): 


kL. kL 
i ke cos] —-sin(@1 } - cos AE 
eae cron] - sof] 


In-e 


2-nurff : 2 
1- sin(61) 
kL. k-L 
ned cos|——-sin]{ 62 ] | — cos} — 
Ine? wie 2 ( ) | 


ony l- sin(62 } : 


E0h(@1) “= j -qo: 


E@hl, "=j “no 


Total E-Field Eh; (01,91): 


E@ht( 61 , $1) := E0h(01 )-AFh(61,61)-4¢ Keheeos(81)_ yc gy yer) Khost 1)} 


E@htl | ‘= [em arn, ¢ -K-h- cos (2) a Ive sad 
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(radians) 


(dimensionless) 


(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


Normalized E-Field E@h;, (81,91): 


E6htn( 61,61) __ [Eeht(61,41)] (V/m) 


max(E@ht1 ) 


Horizontal plane wave reflection coefficient I} (@1) 


: 2 
Th(@1) := costal ea apa sin OL) (dimensionless) 
cos(01) + aferp - sin(@1)° 


aaj __008(02) - jJap- sin(02,)” 


cos (02) 4 ap sin(o2,)” 


¢ component of E-Field of finite dipole Egh(@1): 


(dimensionless) 


gu hat 2 
Eoh(@1) = we - cos KE) (V/m) 
2-nrff 2 
eee . 
Eghl_ :2) wtf - co) (V/m) 
: 2-nrff 2 


Total E-Field Edh; (81,61) & Normalized E-Field Ebh;, (81,61): 


Eght(@1, $1) := Eph(@1 )-AFh(@1 , 61): (<i Kehveos 81) 4 phegy yg) Khost 81)) (V/m) 
j k-h-cos (62 

Eghtl == [evn arn, (e Bren?) | (V/m) 

Epnin(@1 41) := Ean OL. a1)! (v/m) 


max(E¢ht1 ) 
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Radiation Intensity Uh(@1,$61): 
2 
Uh(01 61) sc |Eont(o1 41) 7+ ({E¢htcot, 61) 1)? 
‘yo 


unt [ (leat). (fet 


Radiated Power PV rad: 
2-1 Z 
2 
Phrad .= Uh(61, $1 )-sin( 61) dO1 dol 
0 0 
Phrad = 1.059 
Directivity Dh): 
Uhmax = Uh(¢hp,0) (W/solid angle) Sho ee 
Phrad 
Dho = 5.051 
Effective Isotropic Radiated Power EIRPh: 
EIRPh := Pvrad: Dvo (W) EIRPh = 2.133 
Radiation Resistance RV pt 
Rhr .= —— (Q) Rhr = 135.58 
([In]) 
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(W/solid angle) 


(W/solid angle) 


(dimensionless) 


(dimensionless) 


Input Resistance Rhin: 
Rhin Ses (Q) 


. fi) 
2 


Voltage Relection Coeffecient TH: 


His Rhin- Zo 
Rhin+ Zo 


(dimensionless) 


Reflection Efficiency eve: 


cht :=00d[ 1 - (|rH| | (dimensionless) 


Gain Gh: 
Gh = cht-Dho (dimensionless) 
Gh = 4.633 (dimensionless) 


Polarization Loss Factor PLEh: 


Ehmax := E¢htn(@hp,0) (V/m) 


0 
0 
cha := (dimensionless) 
Ehmax 


lobesh 


PLFh = ets ens Ni (dimensionless) 
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Rhin = 135.58 (Q) 
TH =0.288 (dimensionless) 
cht = 0.917 (dimensionless) 


GhdB =10-log(Gh) (dB) 


GhdB = 6.659 (aB) 


lobesh := floor (+ + } 
X 


PLFh = (dimensionless) 


Maximum Effective Aperture Ahem: 
Gh.w 


4-n 


Ahem ‘= 


-PLFh (m2) Ahem = 35.536 (m2 ) 


Maximum Effective Height nhen: 


hhem := 2, cate (m) hhem =7.15 (m) 
no 
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E-PLANE RADIATION PATTERN: This far-field radiation pattern in the 
(x-z) plane is for the horizontal caged dipole over earth: 


Number of Dipole Conductors: N =8 (elements) 
Radius of Antenna: a=0.125 (meters) 
Height of Antenna Above Earth: h=0.25 (meters) 
: m1 1 : 1 ; 1 
X3(01) := [E@htn} 61 ,—] |-cos} 61+ — Y3(01) ‘= |E@htn] 61 ,—] |-sin} 61 4+ — 
2 Z 2 2 
@ = 0° 


yref 


Y3( 81) 


xref X3¢ 61) 
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H-PLANE RADIATION PATTERN: This far-field radiation pattern in the 
(y-z}) plane is for the horizontal caged dipole over earth: 


Number of Dipole Conductors: N=8 (elements) 
Radius of Antenna: a =0.125 (meters) 
Height of Antenna Above Earth: h=0.25 (meters) 
X2(61) := |Eghtn(61 ,0)| cos + ; Y2(61) := |Eghtn(61 ,0)| sin +3| 
2 2 
6 = 0° 


-- 0.5 
Y2(81) 


xref» X2( 61) 
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THE HELICAL ANTENNA (REFLECTOR OPTION) 
MATHCAD SOFTWARE-DISH_ HEL.MCD 

When built to the proper specifications, the helical antenna possesses 
many qualities which make it suitable for a wide variety of 
communications applications. If the following conditions are satisfied 
the helix will exhibit a highly directional axial main lobe, low side 
lobe level, negligible mutual interference with adjacent antennas, low 
voltage standing wave ratio (VSWR), and resistive input impedance over a 
wide frequency band: 


8 < Cy < 1.15 
n> 3 
12 < a< 14 


(Note: 4 in a subscript indicates the dimension is in wavelengths. 
Mathcad equations can not use symbolic subscripts. Therefore, the symbol 
4} will immediately follow the parameter in equations (i.e., Cd) in lieu 

of subscripts.) 


The helical antenna Mathcad application will compute the following 
parameters (Items with * indicate parameters that are calculated for 
both axial and peripheral feed geometries): 


C = Circumference of Helix 
i = Wavelength 
= Pitch Angle 
Do = Directivity 
p = Relative Phase Velocity 
EQ 4 = Electric Field Components 
Y = Array Factor Phase Shift 
U = Radiation Intensity 
Prad = Radiated Power 
R = Antenna Input Resistance* 
[. = Voltage Reflection Coefficient* 
fe = Reflection Efficiency* 
hem = Maximum Effective Height* 
G = Gain* 
EIRP = Effective Isotropic Radiated Power* 
Rem = Maximum Effective Aperture* 
AR = Axial Ratio 
PLF = Polarization Loss Factor 


BW = Bandwidth 
fhigh = Upper Frequency Limit 
fiow = Lower Frequency Limit 


Acceptable Conductor Diameter 
Ex, y,z = Electric Field Cartesian Coordinates 
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Op + bp = Unit Polarization Vector Coordinate Angles 
Oa = Antenna Unit Polarization Vector 


Imin Mininum Distance to the Far-Field 


The following data must be input based on known or estimated data: 


D = Diameter of Helix (Center to Center) 
S = Spacing Between Turns (Center to Center) 
L = Length Along Conductor of One Turn 
n = Number of Turns 
d = Diameter of Helical Conductor 
f = Frequency of Interest 
m = Desired Mode 
i = Number of Increments in Degrees for Far Field Radiation Pattern 
I, = Antenna Feed Current 
Zo = Characteristic Feed Impedance 
Oy, = Incoming Wave Electric Field Unit Vector 
Enter input data here: 
] 
2 
il) 2 (dimensionless) 

ow =] J 

2 

0 

_ 0.1 : 

D := — (meters) L ‘= 0.097 (meters) 

1 

. ; 10 

n:=10 (turns) f = 3-10 (Hz) 
m ‘= | (dimensionless) Z, ‘= 150 (Q) 
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d = 0.005 
x= ] 
y =1 


(meters) 


(meters) 


(meters) 


(meters) 


z := 1000 


(degrees) 


(A) 


(meters) 


Calculate helical geometric parameters and define constants: 


¢ := 2.9979-10° 
Ci=nD 
C =0.1 
Ch =f 

nN 
Ci = 10.007 
wail 

nN 
Li = 9.70679 
ad 2180. 

tt 
ad = 12.95276 


(meters/sec) 


(meters) 


(meters) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(degrees) 


(degrees) 


No = 120-1 (Q) 
Mee (meters/cycle) 
1 = 9.99310 ° (meters/cycle) 
(2 . 
a := atan|— (radians) 
C 
a = 0.22607 (radians) 
ee. ; : 
Si. iz= (dimensionless) 
d 
Si = 2.30161 (dimensionless) 
kiz22 (m71) 
d 
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Calculate helical antenna parameters: 


Define angular offset 0. Minimum Distance to the Far-Field 
from helical axis: Imin: 
t 

r= 162 (meters) 

Q:= 62295 (radians) 
1 r 16,68 (meters) 
2 
_2-(nS) 
p= 0,(2).2% (radians) oe 1 eae 
1 

rmin ‘= max(r) (meters) 

rmin = 10.58741 (meters) 
Relative Phase Velocity p: 

ad Li 
pat (dimensionless) 
SA + m+ |— 
n 
p =2.89616 (dimensionless) 
Array Factor Phase Shift y. 
: Li ; 
y(8) := 2-n1Si-cos(@) — — (radians) 
p 
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Electric Field Components Eg, Eg: 


E6(6) := E(6) 


E9(6) :=j -E(@) 


Radiation Intensity U(@): 


U(0) = —-({E(0)1)” 
Ty 


Radiated Power Prad: 


2n[2 


U(6)-sin(@) d@ do 


Prad ‘= | 


“0 “0 


Prad = 1.67351°10> 


Directivity Do: 


Do ‘= 12-Ch?-n- Sh (dimensionless) 


Do =2.7658 10° (dimensionless) 


Do2 
Prad 


Do2 = 19.91822 
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._ 47U(0) 


(V/m) 


(V/m) 


(V/m) 


(W/solid angle) 


(dimensionless) 


(dimensionless) 


Axial Ratio AR: 


AR:= ha fine - +)| 
p 


AR = 1.17586 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad-Do (W) 
EIRP = 46.286 (W) 


Polarization Loss Factor PLF: 


| 2 2 
8p -au 7 (radians) 


Zz 


Op = 1.41421 10° (radians) 


Ex =E@(6p)-cos(@p)-cos( 4p) - E¢( 8p )-sin( 4p) 


Ey =E@(8p)-cos( 8p )-sin(p) + Ed @p )-cos( dp) 


Ez =E0(0p)-sin(@p)-- | 


Ex 
a ae 


de [Ex] )?+ (IEyl)?4 (Ezl)? {Ez 


0.5- 0.5} 
oa =| 0.5 + 0.5) 


-9.99999"10 7 


EIRP2 ‘= Prad-Do2 


EIRP2 = 0.03333 


op = anf] 
x 


op =0.7854 
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(dimensionless) 


(dimensionless) 


(radians) 


(radians) 


(V/m) 


(V/m) 


(V/m) 


(dimensionless) 


(dimensionless) 


Polarization Loss Factor PLF: 


PLF = ev ea ) 


(dimensionless) 


Radiation Resistance Ree, 


Prad 
(|Io] )? 


Rr :=2- (Q) 


Dual Parameters 


Axial Feed 


Input Resistance R: 


Ra!=140-4C) (Q) 


Ra = 4.42874 107 (Q) 


Voltage Reflection Coeffecient]: 


Ra- Zo 

Ta -= (dimensionless) 
Ra+Z 
0 

Ta = 0.49399 (dimensionless) 


Reflection Efficiency ¢,: 


era:=1- ([Tal (dimensionless) 


era = 0.75597 (dimensionless) 


174 


PLE = 1 


(dimensionless) 


Rr = 3.34702 10 > (Q) 


Peripheral Feed 


Rp lod: 
dca 
Rp =47.41756 
Rp- Z 
Tp ‘= : 
Rp+ Z, 
Tp =-0.51962 
2 
ep =1-((Tpl) 
erp = 0.72999 


(Q) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain G: 


Ga ‘= era Do 


Gadb := 10-log(era: Do) 


Ga = 2.09088" 10° 


Gadb = 43.20328 


Ga2 := era Do2 


Gadb?2 ‘= 10-log(era:Do2) 


Ga2 = 15.05765 


Gadb2 = 11.77757 


(dimensionless) 


(aB) 


(dimensionless) 


(caB) 


(dimensionless) 


(AaB) 


(dimensionless) 


(aB) 


Maximum Effective Aperture Rem: 


_. fa: i°-Do 
4-n 


Aema PLF 


Aema =0.16615 


2 rah’ Do2- 
4-n 


Aema?2 : PLF 


Aema2 = 1.19657 10° 
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Gp ‘= erp: Do 


Gpdb := 10-log(erp:Do) 


Gp =2.01902-104 


Gpdb = 43.0514 


Gp2 := erp-Do2 


(dimensionless) 


(dB) 


(dimensionless) 


(dB) 


(dimensionless) 


Gpdb2 := 10-log(arp:-Do2) (dB) 


Gp2 = 14.54017 


Gpdb2 = 11.62569 


7 erp-4°-Do. 
4-n 


Aemp: PLF 


Aemp = 0.16044 


ap-h-Do2. 
4-n 


Aemp2 ‘= PLF 


Aemp2 = 1.15545 10 4 


(dimensionless) 


(dB) 


Maximum Effective Height hen: 


hema := 2- JRr- pene (m) hemp := 2- {Rr- Some 

Th No 
hema = 2.42912°10 ° (m) hemp = 2.38701°10 ° 
hema? ‘= 2- |Rr- Bemee (m) hemp2 :=2- {Rr- Bene 

No No 
hema2 =6.51873 10 >” (m) hemp2 =6.40574 10° 
Bandwidth: 
triphone (Hz) flawless 
C 

thigh = 3.44759 10° (Hz) flow = 2.39832: 10° 
BW := fhigh - flow (Hz) 
BW = 1.04927 10° (Hz) 


Acceptable Conductor Diameter: 


dmin ‘= .005-1 (m) dmax := .05-1 


dmin = 4.9965 10° (m) dmax = 4.9965 10 7 
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(Hz) 


(Hz) 


(m) 


HELICAL ANTENNA FAR-FIELD RADIATION PATTERN 


For the purpose of this far-field radiation pattern, the helical antenna 
axis iS equivalent to the Ex = 0 grid line. The pattern is essentially 
symetric when rotated about the antenna's axis. 


Ex(@) := Efe 3] -cos( 9) Ey(6) := efe- 2) si 


0.98 
0.76 


0.54 


Ey(®) 0.1 
0.12 
0.34 
—0.56 


0.78 


PiTAtT TTA TT 


7 
- 
- 
a 
7 
7 
7 
A 
7 


—0.8 0.64 —0.48 —0.32 -0.16 5.551 10 


Ex( 6) 
— EFIELD 
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PARABOLIC REFLECTOR ANALYSIS WITH HELICAL ANTENNA FEED 


The helical antenna Mathcad application (Reflector Option) will compute 
the following additional parameters (Note: a 2 or 22 will follow most 

secondary field parameters in this application, and inorder to speed up 
computation time corresponding vectors are created with subscript A & B: 


Px,y = E-Field Integral and Summation Equation Vector Components 
Fax, ay = Aperture E-Field Components 

@',p = E-Field Integral and Summation Equation Parameters 
Eg = Secondary Field Theta E-Field Component 
Ey = Secondary Field Phi E-Field Component 

U = Secondary Field Radiation Intensity 

Prag = Secondary Field Radiated Power 

Do = Secondary Field Directivity 

G = Secondary Field Gain 

fap = Aperture Efficiency 

fs = Spillover Efficiency 

fe = Taper Efficiency 

i = Phase Efficiency 

€f = Polarization Efficiency 

Eohmic = Ohmic Efficiency 

Ae = Reflector Effective Area 

he = Reflector Effective Height 

EIRP = Reflector Effective Radiated Power 

Ex, y,Z = Secondary Field Cartesian Components 
0, = Reflector Unit Polarization Vector 

PLF = Reflector Polarization Loss Factor 

Imin = Minimum Distance to Secondary Far-Field 
84 = Reflector Periphery Coaltitude 

B+ bp = Theta and Phi for Polarization Loss 


The following additional data must be inputted based on known or 
estimated data: 


E, = Electric Field Scale Factor 

aref = Radius of the Mouth of the Reflector 

foc = Reflector Focal Length 

f) = Blockage Efficiency 

Esp = Spar Efficiency 

tl,t2 = Radiated Power Increments for Phi and Theta 
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i1 = Secondary Field Increments 


N,M = Summation Increments 

repo = Secondary Field Observation Distance 

X,Y,Z = Polarization Loss Factor Coordinates 
Ow = %Incoming Wave Unit Polarization Vector 


Enter additional input data here: 


foc := 5 (meters) 


aref := 5 


(Note: (f£) should be greater 
than minimum distance to the 
far-field for the feed antenna) 


esp = .95 (dimensionless) 
t2 = 10 (increments) 

N =10 (increments) 
rff2 =10° (meters) 


ow2 ‘=| | (dimensionless) 
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eb = .96 
tl =10 
il =10 
M = 10 
x = 1000 
y =1000 
z =1000 


(meters) 


(dimensionless) 


(increments) 


(increments) 


(increments) 


(m) 


Establish Integral and Summation Increments: 


(radians) 


(radians) 


6 


92 --{4)- i ee se eee 
36 il-18 36 36 


Ip = om ace 
t2 


A 
rp], ‘= aref— 
a N 


Calculate reflector geometric parameters: 


00 = 2 atn| Se] (radians) 
2:foc, 
G0 = 0.9273 (radians) 
2 2 
4. 
Sipe ee (m) 
4 foc 
ee rp 
Op) = 2 an (radians) 
2: foc 
4-foc? + (rp1,)? 
pl, :=§———__————_ (m) 
4 foc 
6pl, :=2-at (a (radi ) 
‘= 2-atan radians 
ze 2-f 


A:#1.N 
B:=1.M 
Hl, #20 
#22, = 2-0 
022, {2 AL 4 
13 N 36 
022, := - 
Ty = 1.62 
ro 10-arefi 
ris 8-aref- 
a 


rmin ‘= max(r) 


rmin = 2.0014 10° 
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(increments) 


(increments) 


(radians) 


(radians) 


(radians) 


(radians) 


Aperture Electric Fields Fax: Bayi 


Eax(rp, $1) := (E6( 8p(1p))-- cos( 61) + Ed(@p( rp) )-sin( 41 ))-e! K4 PUP) = foe?) 
Eay(rp, $1) = (E0(@p(rp)):sin( $1) + E6( @p(rp))-cos( 1 ))-e! RCPUP)= fee?) 
Bak! gp :=[E9(p1,,} -(-cos(411,}} + E4(0p1,) -sin(411,) | (pt, — food) 
Eayl, y= (E0(0p1,) -sin(411,) + E¢{6p1,) -cos(411,)) iB OtA~ toed 


Electric Field Integral Equation Vector Components Px, ye 


aref 
Ax(02,62,41) -| Eax(1p,$1)-e sRSepsi 22 oom gh COR) HPL SINE) arg 
0 
2-1 
m2.) | Ax( 62, 62,61) dol 
0 
aref ; 
Ay( 62,42,61 ) -| Eay(rp,$1 )-e -k-rp-sin( 82 )-( cos( 61 )-cos( 62) + sin( >1 )-sin( 92) drp 
0 
2:1 
Py(€2, $2) :| Ay(82, 62,61 ) dol 
0 
ce (223) cost cx(o2e) +sn(1e) sin(29) 
j ‘k-rpl-sin(922 4) - llp)- 22 in{o11p) -sin{ 922, 
Axl, y = Eoxl, ye rpl-sin cos cosl +sin{o sin{@ epi dep 
‘0 
2:4 
Pxl, 3 = Ax], , 461] 
“0 
‘i (23) (oof) cosa) + snr) sinGea) 
j -k-rpl-sin{@22,) - Ip) - 22 in[$l1lp) -sin(o22 
Ayl, ? Eayl, ae rpl-sin cos{ol cosl6 +sin[d sinld rp drp| 
“0 
2-4 
Py aaioe Ay], , d$11 
!0 
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(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


Secondary Field Electric Field Components Eo. Eqi 


i kel em 
£92(62,42) LK (px¢o2,42)-cos(42) + Py(02,62)-sin(42)) (V/m) 
2:nrft2 
j ke KM? 
E§2(02,62) 22S". Px 62, 42)-- sin($2) + Py(02, 42)-cos( $2) )-cos( 62) (V/m) 
2-n-rft2 
joked em 
E03, » eee (Pxl, 300s{422p} +Pyl, a sin(§22,) | (V/m) 
? an > ? 
jeked et? 
E63, 3 = — [Pxl, art sin(§22,)} + Py], -cos{422,) } cos{022,) (V/m) 
UT. 


Secondary Field Radiation Intensity U: 


2 


U2(02,42) = [ le02c02,42)1)"4 (£42062, 42)1)] (W/solid ang) 
My 
U3 .8 = [((e03, 5) (E43, al)’] (W/solid ang) 
1) 


Secondary Field Radiated Power Prag: 
aie n 
Prad2 := (23 4 2 (2.02(025 -sin(622 4 (W) 


Prad2 =0.19318 (W) 


Secondary Field Directivity D,: 


Do22 v= Lt max(U3) (dimensionless) 
Prad2 


Do22 = 1.04762 10° (dimensionless) 
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Spillover Efficiency ¢.>: 


Go 
U(@p)-sin( @p ) dép 
ee 
Es = 
n 
U(6p)-sin( 8p) d@p 
0 
es = 0.19013 


Taper Efficiency ef: 


80 2 


2 
e+ an dip 2eai{ (4) 
Prad 2 2 


et “= ZO 
U(@p)-4-1 
Prad eee 


“0 
et = 0.50837 


Phase Efficiency Epi 


-4.7-6-10° Patee 


ep ‘=e A 


ep = 0.68574 


Aperture Efficiency fap: 


ex [= 98 (dimensionless) eohmic ‘= 98 


cap ‘= eS: et-ep:ex-eb: esp-eohmic (dimensionless) cap =0.05805 
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(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP22 := Prad2:Do22 (Ww) 

EIRP22 = 2.02379 10° (Ww) 

Gain G: 

G22 :=eap-Do22 (dimensionless) G22dB := 10-log(G22) (dimensionless) 
G22 =6.08192 10° (dB) G22dB = 37.84041 (dB) 


Effective Aperture A,: 


Ae22 := eap- 1: aref? (m? ) 


Ae22 = 4.55961 (m2) 


Radiation Resistance Rr: 


Rr22 =2: a (Q) 
([Io] ) 
Rr22 = 0.38636 (Q) 
Effective Height he: 
he22 =2. [Rr22 Ae22 (m) 
No 
he22 =0.13672 (m) 
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Polarization Loss Factor PLF: 


2 2 : 
Op2 and Oe | veacseus! 
Zz 
@p2 = 0.95532 (radians) 
E04 := E62(6p2, p2) (V/m) 


Ex2 ‘= E64-cos( @p2)-cos($p2) - E¢4-sin(¢p2) 


Ey2 ‘= E64-cos(@p2)-sin(¢p2) - Ed4-cos($p2) 


Ez2 


- E64-sin( @p2) 


Ex2 
oa2 ‘= a -|Ey2 


de |Ex2| )? + ([Ey2l )°+ ([Ez2l)? bE22 


0.12034 + 0.48533} 
oa2 =| 0.12034 + 0.48533; 
0.60556 - 0.36501} 


PLF2 := (|owa-saal) 


PLF2 = 0.25003 
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op2 satan (radians) 

op2 =0.7854 (radians) 

E94 ‘= E62(@p2, op2) (V/m) 
(V/m) 


(V/m) 


(V/m) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


PARABOLIC REFLECTOR WITH HELICAL FEED FAR-FIELD ELEVATION PATTERNS 


E-FIELD RADIATION PATTERN: 


"i118 36 


E@2 ( 4 
2 


92 ‘=. Z\- 1 ee SG a ore 
36 


nt 
-sin{ @2 + — 
2 


X2(82) i= 


E62 ( 3] 
° 


cof + 4 Y2(@2) := 


0.032 


0 028 


0 024 


0.021 


¥2( 62) 0017 


0.014 


0.01 


0.003 


0 
—0.003 —0.002 0.002 -0 001 ~0 001 0 0 001 0001 0.002 0.002 0.005 


— EFIELD 
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H-FIELD RADIATION PATTERN: 


-| \ 
X3(02) := [E42(02.0)| cosf0242] 3002) = [Eo2(02.0)| safe 3) 
4 i) 
2} 2 


BEREAN EREe 
tt yt 


0.035 


| 0 028 


0 024 


0 021 


| ¥3( 2) 0.017 


0.014 


0 
0.003 —0 002 “0002 0001 -0 001 0 0.00) 0.001 0 002 0.002 0.005 


“= E FIELD 
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THE SPIRAL ANTENNA (REFLECTOR OPTION 
MATHCAD SOFTWARE-DISH SPI.MCD 


Spiral antennas are a family of two or three dimensional devices that 
possess frequency independant parameters over a wide bandwidth. Spiral 
antennas are commonly used for direction finding, satellite tracking 
and missile guidance. 


The planar spiral may be of the Archemedean, log-spiral, or equiangular 
type. All three radiate two main, circularly polarized lobes 
perpendicular to the plane of the antenna. Additional gain for planar 
spirals may be achieved by placing a metal cavity on the side of the 
antenna with the unwanted lobe. The cavity may be empty or be filled 
with electromagnetic energy absorbing material. These applications 
principly examine the equiangular planar spiral and do not account for 
cavity backed effects. 


(Note: Mathcad equations cannot use symbolic subscripts. Therefore, 
symbols like 7# will immediately follow the parameter in equations. ) 


The spiral antenna Mathcad applications will compute the following 
parameters for eguiangular planar Spirals: 


k = Wavenumber 

~ = Wavelength 

Do = Directivity 

Ed = Electric Field Component 

U = Radiation Intensity 

Umax = Maximum Radiation Intensity 
Prag = Radiated Power 

G = Gain 

EIRP = Effective Isotropic Radiated Power 
Rem = Maximum Effective Aperture 
BW = Bandwidth 

Imin = Minimum Distance to Far-Field 
R, = Radiation Resistance 


Nem = Maximum Effective Height 


fhaigh = Upper Operating Frequency 

flow = Lower Operating Frequency 

ry, = Any point on the nth edge of a spiral 

fex = Equiangular Planar Spiral Expansion Ratio 
Zl = Planar Spiral Input Impedance 

[T = Voltage Reflection Coefficient 

fry = Reflection Efficiency 
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FLF = Polarization Loss Fectcr: 
Mhigh = Upper Operating Wevelensth 
hlow = Lower Operating Wavelength 

A = Planar Spiral Electric Field Amplitude 


The following data must be input based on known or estimated data: 


M Mode 

N Number of Spiral Arms 
f = Frequency of Interest 
1 

xr 

I 


Number of Increments for Far Field Radiation Patterns 


ff = Distance of Far-Field Calculations 
o = Input Current at Antenna Terminals 
Yo = Spiral Feed Point 
a = Flare Rate 
bdn4+] = Angular Arm Width of nth Spiral Arm 
dy = AZimuth to Compute Expansion Ratio 


R = Overall Radius 
Eg = Source Strength Constant for Planar Spirals 


= Wave Unit Polarization Vector 
Zo = Characteristic Impedance of Feed Assembly 


6, = Equiangular Planar Spiral Unit Polarization Vector 


THE PLANAR SPIRAL ANTENNAS 


Enter input data here. 


N22 (arms) f-=3510'° (Hz) 
M ‘= 1 (mode) ; 

lo:=] (amps) 
(Note: Mnayx is N-1) 
1 -= 30 (ancrements) rff-: 110° (meters) 
a‘=.22] (dimensionless) ro = .001 (m) 
Seed wy Eo := 10° (V/m) 
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] 
n 
= 4 
: £2 
6:= (radians) ow c=] (dimensionless) 
. ae 
a 
3a af 
5 Lo | 


Zo ‘= 100 (Q) 


Calculate planar spiral antenna geometric parameters and define 


constants: 


¢ 1229979 108 (meters / sec) Ny = 120-1 (2) 


(meters / cycle) (Farads / m) 


Wy F410 (H / m) 


(dimensionless) 


Calculate planar Spiral antenna parameters : 


Define angular offset 0° 
from y-zZ axis: 


-6n R -6 1 6 . 2a 
10 


Pen (4 pe ee IOS (radians) g:=O— 2a (radians) 
4 , 9 ») i 


“ 1 - 2 
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r 


Distance to Far-Fiele Sas 


Tr) = 16% (m) 
_ &R? 
i —— (m) 
i 
rmin *= max(rr) (m) 
rmin = 0.1 (m) 


Wavenumber Kk: 


7 
kis (m7 2) 
i 


= 


k = 7.33552 10°) (m72) 


Expansion Ratlo fey: 


r(n,ér+ 2:2) 
r(n, 6r) 


eex(n,or) = 


eex( 1,2 2) = 4.00917 


Bandwidth BW: 


Eguiangular Spiral 


Ahigh := 4-ro 
thigh = — 
ahigh 


1, c= 1OR (rm) 


Radial Distance to nth Spiral Edge r: 
(or- Sn ») 


r(n, or) ‘= ro & (m) 


1.22) =40U917* 10> (m) 


(dimensionless) 


(dimensionless) 


Log-Periodic Spiral 


(m) aAhigh| := 20-ro (m) 
Jz: 5 
(Hz) meh high (Hz) 
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thigh =7 49475 107" Ho thigh =~ 49475 10!" Ho, 


Mlow :=4R (m) 


flow =7 49475 10° (Hz) 
BW .= fhigh - flow (Hz) BW] .= thigh] - flow (Hz) 


H 
BW =6 74528 10!° (Hz) BW1 = 74947510" Ae 


Electric Field E(6$) and Electric Field Amplitude A(@): 


wie0 1 (increments) 


Eo-k*-cos(@) (14) a cos(6)) e€ 


diary ae gM 3 Stoo 
' nl ; 


E(6) : — (V/m) 


sin(0)° ff 


9 M (“) atan( a cos(0)) 
cos(a) an e'? 
a 2 
ae) (V/m) 


sin(6) fl + a’ cos( 6)" 


y(8) =v in( 1 +a egeay) + atan( a-cos(6)) (dimensionless) 


fee) 


y- 


E4(0) = AC) MS 


us : 
M (+2) - ao (V/m) 
rff 
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Se seuioral Pe fr en ee tes 
iy ies 2 : & atan{ocosl oo 4) 
2 


n 
col ti (V/m) 
Al = a 
n w Ww : 
snf-3 + 7} 1+ # cof. 2 + | 
2 1 2 1 
Radiation Intensity U(0): 
2A ( 2) 
U(6) = ——-4A(6) (W / solid angle) 
2-1 
Ul = ((a1,)7| (W / solid angle) 
w 2-Ny 
Umax ‘= max( U1) (W / solid angle) 
Umax = 2.2285°10 ° (Ww / solid angle) 


Radiated Power Prag: 


a 
2 

Prad ‘= 4-1: U(6)-sin(6) d6 (W) 
0 

Prad = 7.94528" 10 > (W) 


Directivity Do: 


0 - 4-1 Umax (dimensionless) 
Prad 
Dy = 3.52462 (dimensionless) 
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Effective Isotropic Radiated Power (EIRP): 


cree eS aw) 
= 
EIRP = 0.028 Cw) 


Polarization Loss Factor (PLF): 


—|\2 
PLE = ( low cal) (dimensionless) 


PLF = | (dimensionless) 
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Radiaticn Resistance Fr: Indut Impedance Zi 
4 
ie (Q) — N30a 
2 Zi = ———__— (02) 
(|Io] ) M 
sinfa— 
N 
Rr = 0.01589 (Q) Zi = 1 88496" 10° (Q) 
Voltage Reflection Coefficient [: Reflecton Efficiency €,.: 
pee (dimensionless) om =1-( rl) (dimensionless) 
Zi+ Zo 
T =0.30675 (dimensionless) an =0 9059] (dzmensionless) 
Gain G: 
G:sav D, (dimensionless) GdB ‘= 10 Jog( G) (dB) 
G =3 19298 (dimensionless) GdB = 5 04196 (dB) 


Mazimum Effective Aperture aA 


i” Dy av PLE ; 

cone (m=) 
4a 

Aem = 1 86416 10> (m2 ) 


Maximum Effective Height (hep): 


hem := Rr Aem 2? (m) 
No 
hem = 5 6063°10 ° (m) 
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THE EQUIANGULAR PLANAR SPIRAL ANTENNA FAR-FIELD ELEVATION PATTERN 


For the purpose of this far-field radiation pattern, the spiral antenna 
lies parallel to the Ey = 0 grid line and is centered at the origin. 
The magitude of the electric field pattern 1s rotationally symmetric 
with respect to the Ex=0 grid line. The equiangular planar spiral 
antenna possesses a mirror image radiaton pattern in the -y half plane. 


Egx(6) := |E4(0)| con 3) E¢v(0) := |E4(8)| sof 3 


0.00) 
0.001 


Eoy( 8) 0001 


-4"10 > 


— _ EFIELD 
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PARABOLIC REFLECTOR ANALYSIS WITH SPIRAL ANTENNA FEED 


The spiral antenna Mathcad application (Reflector Option) will compute 
the following additional parameters (Note: a 2 or 22 will follow most 
secondary field parameters in this application, and inorder to speed up 
computation time corresponding vectors are created with subscript A & B: 


Px, y = E-Field Integral and Summation Equation Vector Components 
Fax, ay = Aperture E-Field Components 

@',p = E-Field Integral and Summation Equation Parameters 
Eg = Secondary Field Theta E-Field Component 
E4 = Secondary Field Phi E-Field Component 

U = Secondary Field Radiation Intensity 

Prag = Secondary Field Radiated Power 

De = Secondary Field Directivity 

G = Secondary Field Gain 

fap = Aperture Efficiency 

ft. = Spillover Efficiency 

f& j= Taper Efficiency 

fp = Phase Efficiency 

tf = Polarization Efficiency 

fohmic = Ohmic Efficiency 

Ae = Reflector Effective Area 

he -= Reflector Eitective Height 

EIRP = Reflector Effective Radiated Power 
Ex,y,2z = Secondary Field Cartesian Components 
6, = Reflector Unit Polarization Vector 

PLF = Reflector Polarization Loss Factor 

Imin = Minimum Distance to Secondary Far-Field 
@, = Reflector Periphery Coaltitude 

851 bp = Theta and Phi for Polarization Loss 


The following additional data must be inputted based on known or 
estimated data: 


aref = Radius of the Mouth of the Reflector 

foc = Reflector Focal Length 

ft = Blockage Efficiency 

Esp = Spar Efficiency 

tl,t2 = Radiated Power Increments for Phi and Theta 
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il = Secondary Field Increments 


aref ‘= 0.32 


foc = 0.147 


Summation Increments 
Secondary Field Observation Distance 


Polarization Loss Factor Coordinates 


= Incoming Wave Unit Polarization Vector 


Enter additional input data here: 


(meters) 
(meters) eb = 96 
(dimensionless) tl:=10 
(ancrements) 11 :=10 
(increments) M := 10 
(meters) x:= 1000 
y := 1000 
(dimensionless) 
z := 1000 


198 


(dimensionless) 


(increments) 


(increments) 


(increments) 


(m) 


(m) 


(m) 


Establish Integral and Summation Increments: 


$1 0,07. 24 (radians) A:=1..N (1n¢crements) 
t 
ses Det 3 3 ue 
$2 = 10° ——+ 10°. 2-14 10 (radians) B:=1.M (ancrements) 
1] 
92 :=-(4]- 103, Pe ees $11, 22-02 (radians) 
36 1118 36 36 M 
Ip = 0, HEF ares (m) $22, ‘= 2-ne. (radians) 
t2 
tpl, = are (m) 622, a eae (radians) 
N 18) N 36 
022, :=-— (radians) 
36 
Calculate reflector geometric parameters: 
00 ‘= aan = (radians) Ty: 1.6-2. (m) 
2-foc 
00 = 1.65544 (radians) r<8 10- arefi (m) 
2 2 
orp) 2 4 foc” + rp on a 8-aref (m) 
4 foc ry 
Op(rp) = Data ?P (radians) rmin ‘= max(r) (m) 
2-foc 
2 2 
4-foc’ + [rp] 
ee (m) - rmin = 95.64028 (m) 
4-foc 
1, 
pl, = 2-atan (radians) 
2: foc 
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Aperture Electric Pields Fay, Eay? 
Eax(rp, 61) := Eg @p(rp))-sin( 4] )-e “(PUP fe?) 
Eay(rp, $1 ) "= E¢( @p(rp) )- cos( $1 )-e -k-( pl rp ) = foc: 2) 
) ‘(pla foc: 2) 


Eaxl, p -=E¢(0p! ,)-sin(61 1s): 


Eayl, B= E4(@p! ,)-cos(¢1 I,)-¢ k (p1a~ foe 2) 


Electric Field Integral Equation Vector Components Px, ye 


aref 
Ax(82,62,61) -| Eax(rp,61 )-e k-rp- sim 82 )-( cos( O1)-c0s( $2)-+ sin 91)-8im 92)) 2, ae, 
0 
2a 
Px( 2, 62 ) :| Ax(62, 62,91) dol 
0 
aref . . . 
Ay(02,42,61 ) | Eay(rp, 41 )-e) -k-rp-sin( 82) (cos( $1 )-cos( $2) + sim $1 )-sin( *)) drp 
0 
2-1 
Py( 62, $2) :| Ay(82, $2, 61 ) dol 
0 
= (02) (og) og) sng) (2) 
: j -K-rp]-sin(@22,4)-(cos($1 1p) -cos($22p) + sin (6) 1p) -sin[¢22p 
osl.g| Eaxl, 3 ¢ ‘Tpl 
0 
2-1 
Pxl, 3 AX], dll 
0 
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drp | 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


Secondary Field Electric Field Components Eq. 


E02(62, 42) = SS —_—.(Px(62,42)-cos(42) 4 Py( 02,42 )-sin(62)) 
2 


Pe “kK rff2 


E2( 62, $2) = 1 SS -(Px( 02,42) sin(42) 4 Py( 82 42)-c05(42))-c05(62) 


1 ~ 


eed “ef2 


E03, , a ai cee (Px1, ,-cos(¢22,,) + Pyl, ,-sin(22,)) 


2-n-rff2 


oiB 2-n-rff2 


Secondary Field Radiation Intensity U: 


2 
U2(62, 42) =H [ (E020 02,42)])74 (1E42062.62)1)?] 


2:N 


E83, |)’ + (E03, al)”] 


Prad2 = 1.84163 10 |! 


Secondary Field Directivit Do: 


e 4-1-max(U3) 
Prad2 


Do22: 


Do22 = 4.59318 10° 


[PST ces ( sin($22,)) + Pyl a.p'005($22,) | cos{ 822, } 


201 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(W/solid ang) 


(W/solid ang) 


(W) 


(W) 


(dimensionless) 


(dimensionless) 


Spillover Efficiency co: 


6o 
U(@p)-sin( 8p ) dp 
Hig® 


a Noa 


2 
U(@p)-sin(@p ) d6p 


io 


es = 1.00042 


Taper Efficiency €: 


2 


Bo 


et -= 


et = 0 45135 


Phase Efficiency Epi 


ep = 0.97222 


Aperture Efficiency fap. 


€X -= .98 (dimensionless) eohmic -= .98 


gap -= eS et-ep-ex eb esp cohmic (dimensionless) cap =0.38452 
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(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP22 := Prad2-Do22 


EIRP22 = 8.45893°10 ° 


Gain G: 
G22 °= cap Do22 (dimensionless) G22dB ‘= 10:log( G22) 
G22 =176615*10° (dB) G22dB = 52 47028 


Effective Aperture A,g: 


Ae22 *= eap-n-aref- 


Ae22 = 0.1237 


Radiation Resistance R;: 


Prad2 
([Tol )? 


Rr22 :=2- 


Rr22 = 3.68326 10 |! 


Effective Height he: 


Ae22 
No 


he22 :=2- |Rr22- 


he22 =2.19868 10” 


203 


(W) 


(dimensionless) 


(dB) 


(m2) 


(m2) 


(Q) 


(Q) 


(m) 


(m) 


Polarization Loss Factor PLF: 


[2,2 | 
@p2 := an( 227 (radians) 


z 


@p2 = 0.95532 (radians) 
E64 := E62(@p2, ¢p2) (V/m) 

Ex2 := E64-cos(@p2)-cos( 6p2) — E$4-sin(¢p2) 
Ey2 ‘= E04-cos(@p2)-sin(¢p2) - E$4-cos($p2) 


Ez2 := - E64 sin( 6p2) 


x2 
oa2 ie ee Le ee Ey2 


dc tex2t y+ CIEy2I 74 ({Ez2l)? \Ez2 


0.40672 + 0.03531) 
oa2 = |-0.40672 + 0.03531) 
0.81344 - 0.07063; 


PLF2 := (|ow2'saal) 


PLF2 = 0.16667 
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$p2 san(2| (radians) 

¢p2 = 0.7854 (radians) 

E64 := Eo2(@p2, ¢p2) (V/m) 
(V/m) 


(V/m) 


(V/m) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


PARABOLIC REFLECTOR WITH SPIRAL FEED ANTENNA 


E-PANE RADIATION PATTERN: 


il :=20 an (2]- [Oe ise ee Oe oa 
1118 36 36 


X1(02) *= ee2(oe. 3] 


cos{ea :| Y1(02) ‘= ev2(ee.2] 55 


3.09 10 15 


3.08910 !9 


3.08710 !° 


3.086 10 !9 


3.08410 15 


Y1( 62) 3.083 10 !9 


“OA 
YE ET 


-15 


3.08 10 


3.07% 10 


3.077 10 !9 


3.076 10 19 


— EFIELD 
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H-PLANE RADIATION PATTERN: 


X2(02) := |E42(62,0)| coo + ;) Y2(02) := {E42(02,0)| snl 2 - 3| 


Tt 
PEt int AT 
SEA 
LUTE TT 


1310 ° 


1.18 10° 


1.06910 


~7 


9 4°10 


8.2 10” 


Y @2) 710 7 


5.8710 ” 


4610? 


— EFIELD 
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10. 


Is Le 


12. 


13. 
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